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Abstract 

Passive linear networks are used in a wide variety of engineering applications, but the best 
studied are electrical circuits made of resistors, inductors and capacitors. We describe a cat¬ 
egory where a morphism is a circuit of this sort with marked input and output terminals. In 
this category, composition describes the process of attaching the outputs of one circuit to the 
inputs of another. We construct a functor, dubbed the ‘black box functor’, that takes a circuit, 
forgets its internal structure, and remembers only its external behavior. Two circuits have the 
same external behavior if and only if they impose same relation between currents and potentials 
at their terminals. The space of these currents and potentials naturally has the structure of a 
symplectic vector space, and the relation imposed by a circuit is a Lagrangian linear relation. 
Thus, the black box functor goes from our category of circuits to the category of symplectic vec¬ 
tor spaces and Lagrangian linear relations. We prove that this functor is a symmetric monoidal 
dagger functor between dagger compact categories. We assume the reader has some familiarity 
with category theory, but none with circuit theory or symplectic linear algebra. 
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1 Introduction 


In late 1940s, just as Feynman was developing his diagrams for processes in particle physics, Eilenberg 
and Mac Lane initiated their work on category theory. Over the subsequent decades, and especially 
in the work of Joyal and Street in the 1980s [25, 26], it became clear that these developments were 
profoundly linked: monoidal categories have a precise graphical representation in terms of string 
diagrams, and conversely monoidal categories provide an algebraic foundation for the intuitions 
behind Feynman diagrams. The key insight is the use of categories where morphisms describe 
physical processes, rather than structure-preserving maps between mathematical objects [7, 14]. 

In work on fundamental physics, the cutting edge has moved from categories to higher cate¬ 
gories [6]. But the same techniques have filtered into more immediate applications, particularly in 
computation and quantum computation [1, 3, 40]. This paper is part of a still nascent program of 
applying string diagrams to engineering, with the aim of giving diverse diagram languages a unified 
foundation based on category theory [4, 10, 27, 37, 43]. 

Indeed, even before physicists began using Feynman diagrams, various branches of engineering 
were using diagrams that in retrospect are closely related. Foremost among these are the ubiquitous 
electrical circuit diagrams. Although less well-known, similar diagrams are used to describe networks 
consisting of mechanical, hydraulic, thermodynamic and chemical systems. Further work, pioneered 
in particular by Forrester [22] and Odum [33], applies similar diagrammatic methods to biology, 
ecology, and economics. 

As discussed in detail by Olsen [34], Paynter [35] and others, there are mathematically precise 
analogies between these different systems. In each case, the system’s state is described by variables 
that come in pairs, with one variable in each pair playing the role of ‘displacement’ and the other 
playing the role of ‘momentum’. In engineering, the time derivatives of these variables are sometimes 
called ‘flow’ and ‘effort’. In classical mechanics, this pairing of variables is well understood using 
symplectic geometry. Thus, any mathematical formulation of the diagrams used to describe networks 
in engineering needs to take symplectic geometry as well as category theory into account. 



displacement 

flow 

momentum 

effort 


q 

q 

P 

P 

Electronics 

charge 

current 

flux linkage 

voltage 

Mechanics (translation) 

position 

velocity 

momentum 

force 

Mechanics (rotation) 
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angular velocity 

angular momentum 
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Hydraulics 

volume 

flow 

pressure momentum 

pressure 

Thermodynamics 

entropy 

entropy flow 

temperature momentum 

temperature 

Chemistry 

moles 

molar flow 

chemical momentum 

chemical potential 


While diagrams of networks have been independently introduced in many disciplines, we do not 
expect formalizing these diagrams to immediately help the practitioners of these disciplines. At 
hrst the flow of information will mainly go in the other direction: by translating ideas from these 
disciplines into the language of modern mathematics, we can provide mathematicians with food for 
thought and interesting new problems to solve. We hope that in the long run mathematicians can 
return the favor by bringing new insights to the table. 

Although we shall keep the broad applicability of network diagrams in the back of our minds, we 
couch our discussion in terms of electrical circuits, for the sake of familiarity. In this paper our goal 
is somewhat limited. We only study circuits built from ‘passive’ components: that is, those that do 
not produce energy. Thus, we exclude batteries and current sources. We only consider components 
that respond linearly to an applied voltage. Thus, we exclude components such as nonlinear resistors 
or diodes. Finally, we only consider components with one input and one output, so that a circuit 
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can be described as a graph with edges labeled by components. Thus, we also exclude transformers. 
The most familiar components our framework covers are linear resistors, capacitors and inductors. 

While we hope to expand our scope in future work, the class of circuits made from these com¬ 
ponents has appealing mathematical properties, and is worthy of deep study. Indeed, this class has 
been studied intensively for many decades by electrical engineers [2, 13, 41]. Even circuits made 
exclusively of resistors have inspired work by mathematicians of the caliber of Weyl [46] and Smale 
[ 42 ]. 

Our work relies on this research. All we are adding is an emphasis on symplectic geometry 
and an explicitly ‘compositional’ framework, which clarifies the way a larger circuit can be built 
from smaller pieces. This is where monoidal categories become important: the main operations for 
building circuits from pieces are composition and tensoring. 

Our strategy is most easily illustrated for circuits made of linear resistors. Such a resistor dissi¬ 
pates power, turning useful energy into heat at a rate determined by the voltage across the resistor. 
However, a remarkable fact is that a circuit made of these resistors always acts to minimize the 
power dissipated this way. This ‘principle of minimum power’ can be seen as the reason symplectic 
geometry becomes important in understanding circuits made of resistors, just as the principle of 
least action leads to the role of symplectic geometry in classical mechanics. 

Here is a circuit made of linear resistors: 



The wiggly lines are resistors, and their resistances are written beside them: for example, 3H means 
3 ohms, an ‘ohm’ being a unit of resistance. To formalize this, define a circuit of linear resistors to 
consist of: 

• a set of nodes, 

• a set A of edges, 

• maps s, t: E ^ N sending each edge to its source and target node, 

• a map r: E ^ (0,oo) specifying the resistance of the resistor labelling each edge, 

• maps i: X ^ N, o: Y N specifying the inputs and outputs of the circuit. 

When we run electric current through such a circuit, each node n G N gets a ‘potential’ (j){n). 
The ‘voltage’ across an edge e G E is defined as the change in potential as we move from to the 
source of e to its target, (t>{t{e)) — ^(s(e)), and the power dissipated by the resistor on this edge 
equals 

The total power dissipated by the circuit is therefore twice 

e^E ^ ^ 

The factor of 4 is convenient in some later calculations. Note that P is a nonnegative quadratic 
form on the vector space However, not every nonnegative definite quadratic form on arises 
in this way from some circuit of linear resistors with N as its set of nodes. The quadratic forms that 
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do arise are called ‘Dirichlet forms’. They have been extensively investigated [23, 30, 38, 39], and 
they play a major role in our work. 

We write ON = i{X) U o{Y) for the set of ‘terminals’: that is, nodes corresponding to inputs 
and outputs. The principle of minimum power says that if we fix the potential at the terminals, the 
circuit will choose the potential at other nodes to minimize the total power dissipated. An element 
ip of the vector space assigns a potential to each terminal. Thus, if we fix ip, the total power 
dissipated will be twice 

Q{ip) = min P{(p) 

The function Q : R is again a Dirichlet form. We call it the ‘power functional’ of the circuit. 

Now, suppose we are unable to see the internal workings of a circuit, and can only observe its 
‘external behavior’: that is, the potentials at its terminals and the currents flowing into or out of 
these terminals. Remarkably, this behavior is completely determined by the power functional Q. 
The reason is that the current at any terminal can be obtained by differentiating Q with respect 
to the potential at this terminal, and relations of this form are all the relations that hold between 
potentials and currents at the terminals. 

The Laplace transform allows us to generalize this immediately to circuits that can also contain 
linear inductors and capacitors, simply by changing the field we work over, replacing R by the field 
R(s) of rational functions of a single real variable, and talking of ‘impedance’ where we previously 
talked of resistance. We obtain a category Circ where, roughly speaking, an object is a finite set, 
a morphism f: X T is a circuit with input set X and output set D, and composition is given 
by identifying the outputs of one circuit with the inputs of the next, and taking the resulting union 
of labelled graphs. Each such circuit gives rise to a Dirichlet form, now defined over R(s), and this 
Dirichlet form completely describes the externally observable behavior of the circuit. 

We can take equivalence classes of circuits, where two circuits count as the same if they have the 
same Dirichlet form. We wish for these equivalence classes of circuits to form a category. Although 
there is a notion of composition for Dirichlet forms, we find that it lacks identity morphisms or, 
equivalently, it lacks morphisms representing ideal wires of zero impedance. To address this we turn 
to Lagrangian subspaces of symplectic vector spaces. These generalize quadratic forms via the map 

(q : ^ f) ^ (Graph(dQ) = {{^p, dQ^) \ ^p € C F®^ © 

taking a quadratic form Q on the vector space F^^ over a field F to the graph of its differential dQ. 
Here we think of the symplectic vector space F^-^ © (F^-^)* as the state space of the circuit, and the 
subspace Graph(dQ) as the subspace of attainable states, with ip G F^^ describing the potentials at 
the terminals, and dQ^ G (F^^)* the currents. 

This construction is well known in classical mechanics [45], where the principle of least action 
plays a role analogous to that of the principle of minimum power here. The set of Lagrangian 
subspaces is actually an algebraic variety, the ‘Lagrangian Grassmannian’, which serves as a com- 
pactification of the space of quadratic forms. The Lagrangian Grassmannian has already played a 
role in Sabot’s work on circuits made of resistors [38, 39]. For us, its importance it that we can 
find identity morphisms for the composition of Dirichlet forms by taking circuits made of parallel 
resistors and letting their resistances tend to zero: the limit is not a Dirichlet form, but it exists 
in the Lagrangian Grassmannian. Indeed, there exists a category LagrRel with finite dimensional 
symplectic vector spaces as objects and ‘Lagrangian relations’ as morphisms: that is, linear relations 
from V to W that are given by Lagrangian subspaces of E © W, where V is the symplectic vector 
space conjugate to V. 

To move from the Lagrangian subspace defined by the graph of the differential of the power 
functional to a morphism in the category LagrRel—that is, to a Lagrangian relation—we must treat 
seriously the input and output functions of the circuit. These express the circuit as built upon a 
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cospan 


N 



X Y. 


Applicable far more broadly than this present formalization of circuits, cospans model systems with 
two ‘ends’, an input and output end, albeit without any connotation of directionality: we might just 
as well exchange the role of the inputs and outputs by taking the mirror image of the above diagram. 
The role of the input and output functions, as we have discussed, is to mark the terminals we may 
glue onto the terminals of another circuit, and the pushout of cospans gives formal precision to this 
gluing construction. 

One upshot of this cospan framework is that we may consider circuits with elements of N that 
are both inputs and outputs, such as this one: 



This corresponds to the identity morphism on the finite set with two elements. Another is that some 
points may be considered an input or output multiple times; we draw this: 


inputs 



<-outputs 


This allows us to connect two distinct outputs to the above double input. 

Given a set X of inputs or outputs, we understand the electrical behavior on this set by con¬ 
sidering the symplectic vector space © (F^) , the direct sum of the space F^ of potentials and 
the space (F^) of currents at these points. A Lagrangian relation specifies which states of the 
output space F^ © (F^) are allowed for each state of the input space F^ © (F^) . Turning the 
Lagrangian subspace Graph((i(5) of a circuit into this information requires that we understand the 
‘symplectification’ 

S'/: F^ © (F®)* ^ F^ © (F^)* 

and ‘twisted symplectification’ 

SV: F® © (F^)* ^ F^ © (F^)* 

of a function f: A ^ B between finite sets. In particular we need to understand how these apply 
to the input and output functions with codomain restricted to dN; abusing notation, we also write 
these z: A —>■ dN and o: Y ^ dN. 

The symplectification is a Lagrangian relation, and the catch phrase is that it ‘copies voltages’ 
and ‘splits currents’. More precisely, for any given potential-current pair in F® © (F-®) , its 

image under Sf comprises all elements of {tp', l') G F^ © (F^) such that the potential at a G A is 
equal to the potential at /(a) G B, and such that, for each fixed b € B, collectively the currents 
at the a G f~^{b) sum to the current at b. We use the symplectification So of the output function 
to relate the state on dN to that on the outputs Y. As our current framework is set up to report 
the current out of each node, to describe input currents we define the twisted symplectification 
S^f: F^ © (F^)* —>• F'^^ © (F^)* almost identically to the above, except that we flip the sign of the 
currents l' G (F"^^)*. We use the twisted symplectification S*i of the input function to relate the 
state on dN to that on the inputs. 
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The Lagrangian relation corresponding to a circuit is then the set of all potential-current pairs 
that are possible at the inputs and outputs of that circuit. For instance, consider a resistor of 
resistance r, with one end considered as an input and the other as an output: 

r 

input-► •-^-• < -output 


To obtain the corresponding Lagrangian relation, we must first specify domain and codomain sym- 
plectic vector spaces. In this case, as the input and output sets each consist of a single point, these 
vector spaces are both F © F*, where the first summand is understood as the space of potentials, 
and the second the space of currents. 

Now, the resistor has power functional Q: F^ —)■ F given by 

Qi'lpl,'lp2) = ^(V’2 - V’l)^ 

and the graph of the differential of Q is 

Graph(dQ) = - ^>2), y(V'2 - V'l)) ki, ^^2 e F} C f 2 © (F 2 )T 

In this example the input and output functions i, o are simply the identity functions on a one element 
set, so the symplectification of the output function is simply the identity linear transformation, and 
the twisted symplectification of the input function is the isomorphism between conjugate symplectic 
vector spaces F © F* —>■ F © F* mapping (</>, i) to {(j), —i) This implies that the behavior associated 
to this circuit is the Lagrangian relation 

{(V'i,b V’2,*) I e F,i = - V'l)} c F ©F* ©F©F*. 

This is precisely the set of potential-current pairs that are allowed at the input and output of 
a resistor of resistance r. In particular, the relation i = ~ ' 4 ’ i) is well known in electrical 

engineering: it is ‘Ohm’s law’. 

A crucial fact is that the process of mapping a circuit to its corresponding Lagrangian relation 
identifies distinct circuits. For example, a single 2-ohm resistor: 

2n 

input- > • -^-• <-output 


has the same Lagrangian relation as two I-ohm resistors in series: 


input- > 




4-output 


The Lagrangian relation does not shed any light on the internal workings of a circuit. Thus, we call 
the process of computing this relation ‘black boxing’: it is like encasing the circuit in an opaque box, 
leaving only its terminals accessible. Fortunately, the Lagrangian relation of a circuit is enough to 
completely characterize its external behavior, including how it interacts when connected with other 
circuits. 

Put more precisely, the black boxing process is functorial: we can compute the black boxed 
version of a circuit made of parts by computing the black boxed versions of the parts and then 
composing them. In fact we shall prove that Circ and LagrRel are dagger compact categories, and 
the black box functor preserves all this structure: 

Theorem 1.1. There exists a symmetric monoidal dagger functor, the black box functor 

■: Circ —> LagrRel, 
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mapping a finite set X to the symplectic vector space © (F^)* it generates, and a circuit ((TV, E, s, t, r), i, o) 
to the Lagrangian relation 

U S\{v) X So{v) C F^ © (F^)* ©F^ © (F^)*, 

'uGGraph(rfQ) 

where Q is the circuit’s power functional. 

The goal of this paper is to prove and explain this result. The proof itself is more tricky than one 
might first expect, but our approach introduces various concepts that should be useful throughout the 
study of networks, such as ‘decorated cospans’ and ‘corelations’. These provide a general framework 
for discussing open networked systems—and not only the passive linear systems discussed here, but 
also others, such as Markov processes [5]. 

Cospans are already familiar as a formalism for making entities with an arbitrarily designated 
‘input end’ and ‘output end’ into the morphisms of a category. For example, in topological quantum 
field theory we use special cospans called ‘cobordisms’ to describe pieces of spacetime [6, 7]. In 
Section 4 we introduce ‘decorated cospans’ to describe circuits with specified inputs and outputs. 

Later, with the help of machinery developed in a companion paper [20], we use decorated cospans 
to set up several functors that appear as factors of the black box functor, as steps in proving its 
functoriality. 

Just as a relation is an isomorphism class of spans of a special sort, namely jointly monic spans, 
a ‘corelation’ is an isomorphism class of jointly epic cospans. While corelations are less widely used 
than relations, they turn out to be perfectly suited for describing circuits of ideal perfectly conductive 
wires—precisely the class of circuits that cannot be modeled simply as Dirichlet forms, since the 
power functional becomes infinite when it includes terms where the resistance is zero. We introduce 
corelations in Section 7, and they too play a key role in constructing the black box functor. 

With these tools in hand, the black box functor turns out to rely on a tight relationship between 
Kirchhoff’s laws, the minimization of Dirichlet forms, and the ‘symplectification’ of corelations. It 
is well known that away from the terminals, a circuit must obey two rules known as Kirchhoff’s 
laws. We have already noted that the principle of minimum power states that a circuit will ‘choose’ 
potentials on its interior that minimize the power functional. We clarify the relation between these 
points in Theorems 3.4 and 3.5, which together show that minimizing a Dirichlet form over some 
subset amounts to assuming that the corresponding circuit obeys Kirchhoff’s laws on that subset. 

We have also mentioned the symplectification of functions above. Extending this to allow sym¬ 
plectification of corelations, this process gives a map sending corelations to Lagrangian relations that 
describe the behavior of ideal perfectly conductive wires. We prove that these symplectified core¬ 
lations simultaneously impose Kirchhoff’s laws (Proposition 7.7 and Example 7.8) and accomplish 
the minimization of Dirichlet forms (Theorem 8.3). 

Together, our results show that these three concepts—Kirchhoff’s laws from circuit theory, 
the analytic idea of minimizing power dissipation, and the algebraic idea of symplectification of 
corelations—are merely different faces of one law: the law of composition of circuits. 

1.1 Finding your way through this paper 

This paper is split into three parts, addressing in turn the questions: 

1. What do circuit diagrams mean? 

11. How do we interact with circuit diagrams? 

III. How is meaning preserved under these interactions? 

We begin Part I, on the semantics of circuit diagrams, with a discussion of circuits of linear 
resistors, developing the intuition for the governing laws of passive linear circuits—Ohm’s law. 



Kirchhoff’s voltage law, and Kirchhoff’s current law—in a time-independent setting (Section 2). 
This allows us to develop the concept of Dirichlet form as a representation of power consumption, 
and understand their composition as minimizing power, an expression of the current law. In Section 
3, the Laplace transform then allows us to recapitulate these ideas after introducing inductors and 
capacitors, speaking of impedance where we formerly spoke of resistance, and generalizing Dirichlet 
forms from the field R to the field M(s) of real rational functions. While in this setting the principle of 
minimum power is replaced by a variational principle, the intuitions gained from circuits of resistors 
still remain useful. 

At the end of this part, we show that Dirichlet forms alone do not provide the flexibility to 
construct a category representing the semantics of circuit diagrams. The fundamental reason is that 
they do not allow us to describe circuits involving ideal perfectly conductive wires. This motivates 
the development of more powerful machinery. 

Part II, on the syntax of circuit diagrams, contains the main technical contributions of the 
paper. It begins with Section 4, which develops machinery to construct what we call categories of 
decorated cospans. These are categories where the objects are finite sets and the morphisms are 
cospans in the category of finite sets together with some extra structure on the apex. Circuits, as 
defined above, are naturally an example of such a construction, and Section 5 lays out the details 
of this. In Section 6 we then review the basic theory of linear Lagrangian relations, giving details 
to the correspondence we have defined between Dirichlet forms, and hence passive linear circuits, 
and Lagrangian relations. Section 7 then takes immediate advantage of the added flexibility of 
Lagrangian relations, discussing the ‘trivial’ circuits comprising only perfectly conductive wires, 
which mediate the notion of composition of circuits. 

Having introduced these prerequisites, we get to the point in Part III. In Section 8 we introduce 
the black box functor, and in Section 9 we prove our main result. 
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Part I 

Passive Linear Circuits 


In this part we review the meaning of circuit diagrams comprising resistors, inductors, and capacitors, 
giving an answer to the question “What do circuit diagrams mean?”. 

To elaborate, while circuit diagrams model electric circuits according to their physical form, 
another, often more relevant, way to understand a circuit is by its external behavior. This means 
the following. To an electric circuit we associate two quantities to each edge: voltage and current. 
We are not free, however, to choose these quantities as we like; circuits are subject to governing laws 
that imply voltages and currents must obey certain relationships. From the perspective of control 
theory we are particularly interested in the values these quantities take at the so-called terminals, 
and how altering one value will affect the other values. We call two circuits equivalent when they 
determine the same relationship. Our main task in this first part is to explore when two circuits are 
equivalent. 
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2 Circuits of linear resistors 


In order to let physical intuition lead the way, we begin by specialising to the case of linear resistors. 
In this section we describe how to find the function of a circuit from its form, advocating in particular 
the perspective of the principle of minimum power. This allows us to identify the external behavior 
of a circuit with a so-called Dirichlet form representing the dependence of its power consumption on 
potentials at its terminals. 

2.1 Circuits as labelled graphs 

The concept of an abstract open electrical circuit made of linear resistors is well known in electrical 
engineering, but we shall need to formalize it with more precision than usual. The basic idea is 
that a circuit of linear resistors is a graph whose edges are labelled by positive real numbers called 
‘resistances’, and whose sets of vertices is equipped with two subsets: the ‘inputs’ and ‘outputs’. 
This unfolds as follows. 

A (closed) circuit of resistors looks like this: 



We can consider this a labelled graph, with each resistor an edge of the graph, its resistance its label, 
and the vertices of the graph the points at which resistors are connected. 

A circuit is ‘open’ if it can be connected to other circuits. To do this we hrst mark points at 
which connections can be made by denoting some vertices as input and output terminals: 



Then, given a second circuit, we may choose a relation between the output set of the first and the 
input set of this second circuit, such as the simple relation of the single output vertex of the circuit 
above with the single input vertex of the circuit below. 



We connect the two circuits by identifying output and input vertices according to this relation, 
giving in this case the composite circuit: 
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More formally, we define a graph'^ to be a pair of functions s,t: E —> N where E and N are 
finite sets. We call elements of E edges and elements of N vertices or nodes. We say that the 
edge e € E has source s(e) and target t(e), and also say that e is an edge from s(e) to t(e). 

To study circuits we need graphs with labelled edges: 

Definition 2.1. Given a set L of labels, an T-graph is a graph equipped with a function r: E ^ L: 

L E N. 

t 

For circuits made of resistors we take L = (0, oo), but later we shall take L to be a set of positive 
elements in some more general field. In either case, a circuit will be an L-graph with some extra 
structure: 

Definition 2.2. Given a set L, a circuit over L is an L-graph L i - E { N together with 

t 

finite sets X, Y, and functions i: X ^ N and o: Y N. We call the sets i{X), o(Y), and 
dN = i{X) U o{Y) the inputs, outputs, and terminals or boundary of the circuit, respectively. 

We will later make use of the notion of connectedness in graphs. Recall that given two vertices 
v,w S of a graph, a path from t; to ic is a finite sequence of vertices v = vo,vi,... ,Vn = w 
and edges ei,..., e„ such that for each 1 < * < n, either is an edge from Vi to or an edge 
from Ui+i to Vi. A subset S of the vertices of a graph is connected if, for each pair of vertices in S, 
there is a path from one to the other. A connected component of a graph is a maximal connected 
subset of its vertices.^ 

In the rest of this section we take L = (0, oo) C R and fix a circuit over (0, oo). The edges of 
this circuit should be thought of as ‘wires’. The label re G (0, oo) stands for the resistance of the 
resistor on the wire e. There will also be a voltage and current on each wire. In this section, these 
will be specified by functions V € R.^ and I € R.^. Here, as customary in engineering, we use I for 
‘intensity of current’, following Ampere. 

2.2 Ohm’s law, Kirchhoff’s laws, and the principle of minimum power 

In 1827, Georg Ohm published a book which included a relation between the voltage and current 
for circuits made of resistors [32]. At the time, the critical reception was harsh: one contemporary 
called Ohm’s work “a web of naked fancies, which can never find the semblance of support from even 
the most superficial of observations”, and the German Minister of Education said that a professor 
who preached such heresies was unworthy to teach science [17, 24]. However, a simplified version of 
his relation is now widely used under the name of ‘Ohm’s law’. We say that Ohm’s law holds if 
for all edges e G E the voltage and current functions of a circuit obey: 

V{e) = r(e)/(e). 

Kirchhoff’s laws date to Gustav Kirchhoff in 1845, generalising Ohm’s work. They were in turn 
generalized into Maxwell’s equations a few decades later. We say Kirchhoff’s voltage law holds 
if there exists (p G such that 

V(e) = cP(t(e)) -<p(s(e)). 

We call the function cj) a potential, and think of it as assigning an electrical potential to each node 
in the circuit. The voltage then arises as the differences in potentials between adjacent nodes. If 

^In this paper we refer to directed multigraphs simply as graphs. 

^In the theory of directed graphs the qualifier ‘weakly’ is commonly used before the word ‘connected’ in these two 
definitions, in distinction from a stronger notion of connectedness requiring paths to respect edge directions. As we 
never consider any other sort of connectedness, we omit this qualifier. 
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Kirchhoff’s voltage law holds for some voltage V, the potential <j) is unique only in the trivial case of 
the empty circuit: when the set of nodes N is empty. Indeed, two potentials define the same voltage 
function if and only if their difference is constant on each connected component of the graph. 

We say Kirchhoflf’s current law holds if for all nonterminal nodes n G N \ dN we have 



s(e)=n t{e)—n 

This is an expression of conservation of charge within the circuit; it says that the total current 
flowing in or out of any nonterminal node is zero. Even when Kirchhoff’s current law is obeyed, 
terminals need not be sites of zero net current; we call the function l G that takes a terminal 
to the difference between the outward and inward flowing currents, 

L-.dN —^ M 

t{e)—n s{e)—n 

the boundary current for I. 

A boundary potential is also a function in but instead thought of as specifying potentials 
on the terminals of a circuit. As we think of our circuits as open circuits, with the terminals points 
of interaction with the external world, we shall think of these potentials as variables that are free for 
us to choose. Using the above three principles—Ohm’s law, Kirchhoff’s voltage law, and Kirchhoff’s 
current law—it is possible to show that choosing a boundary potential determines unique voltage 
and current functions on that circuit. 

The so-called ‘principle of minimum power’ gives some insight into how this occurs, by describing 
a way potentials on the terminals might determine potentials at all nodes. From this, Kirchhoff’s 
voltage law then gives rise to a voltage function on the edges, and Ohm’s law gives us a current 
function too. We shall show, in fact, that a potential satisfies the principle of minimum power for a 
given boundary potential if and only if this current obeys Kirchhoff’s current law. 

A circuit with current / and voltage V dissipates energy at a rate equal to 

^/(e)U(e). 

eeE 

Ohm’s law allows us to rewrite I as V/r, while Kirchhoff’s voltage law gives us a potential (j) such 
that V{e) can be written as (j){t{e)) — (j){s{e)), so for a circuit obeying these two laws the power can 
also be expressed in terms of this potential. We thus arrive at a functional mapping potentials </> to 
the power dissipated by the circuit when Ohm’s law and Kirchhoff’s voltage law are obeyed for (j). 

Definition 2.3. The extended power functional P: > K. 0 / a circuit is defined by 

eeE '' ' 

The factor of ^ is inserted to cancel the factor of 2 that appears when we differentiate this expression. 
We call P the extended power functional as we shall see that it is defined even on potentials that are 
not compatible with the three governing laws of electric circuits. We shall later restrict the domain 
of this functional so that it is defined precisely on those potentials that are compatible with the 
governing laws. Note that this functional does not depend on the directions chosen for the edges of 
the circuit. 

This expression lets us formulate the ‘principle of minimum power’, which gives us information 
about the potential (j) given its restriction to the boundary of T. Call a potential G an 
extension of a boundary potential ip G if is equal to ip when restricted to —that is, if 

(pldN = V’- 
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Definition 2.4. We say a potential (j> G obeys the principle of minimum power for a 

boundary potential if G if (f minimizes the extended power functional P subject to the constraint 
that 4> is an extension of if. 

As promised, in the presence of Ohm’s law and Kirchhoff’s voltage law, the principle of minimum 
power is equivalent to Kirchhoff’s current law. 

Proposition 2.5. Let (f be a potential extending some boundary potential if. Then (f obeys the 
principle of minimum power for if if and only if the current 

/(e) = -^{(f{t{e)) - (f{s{e))) 
r(e) 


obeys Kirchhoff’s current law. 

Proof. Fixing the potentials at the terminals to be those given by the boundary potential if, the 
power is a nonnegative quadratic function of the potentials at the nonterminals. This implies that 
an extension (f of if minimizes P precisely when 

dP{T) 

d(p{n) 

for all nonterminals n G N \ dN. Note that the partial derivative of the power with respect to the 
potential at n is given by 



dP 


dip{n) 


v=4‘ 


Y. -((.(s(e))) 

t{e)—n 

Y. Y 

t{e)—n s{e)—n 


Y We))-((.(s(e))) 

s{e)=n ^ ' 


Thus (f obeys the principle of minimum power for if if and only if 

^ 7(e) = ^ 7(e) 

s(e)=n t(e)=n 


for all n G N \ dN, and so if and only if Kirchhoff’s current law holds. 


□ 


2.3 A Dirichlet problem 

We remind ourselves that we are in the midst of understanding circuits as objects that define 
relationships between boundary potentials and boundary currents. This relationship is defined by 
the stipulation that voltage-current pairs on a circuit must obey Ohm’s law and Kirchhoff’s laws— 
or equivalently. Ohm’s law, Kirchhoff’s voltage law, and the principle of minimum power. In this 
subsection we show these conditions imply that for each boundary potential if on the circuit there 
exists a potential (f on the circuit extending if, unique up to what may be interpreted as a choice of 
reference potential on each connected component of the circuit. From this potential (f we can then 
compute the unique voltage, current, and boundary current functions compatible with the given 
boundary potential. 

Fix again a circuit with extended power functional P: —>■ R. Let V: —>■ be the 

operator that maps a potential </> G R.^ to the function from to R given by 


dP 

dip{n) 


v=4> 
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As we have seen, this function takes potentials to twice the pointwise currents that they induce. We 
have also seen that a potential (j) is compatible with the governing laws of circuits if and only if 


— O' 


( 1 ) 


The operator V acts as a discrete analogue of the Laplacian for the graph T, so we call this operator 
the Laplacian of T, and say that equation (1) is a version of Laplace’s equation. We then say that 
the problem of finding an extension (j) of some fixed boundary potential tp that solves this Laplace’s 
equation—or, equivalently, the problem of finding a (p that obeys the principle of minimum power 
for Ip —is a discrete version of the Dirichlet problem. 

As we shall see, this version of the Dirichlet problem always has a solution. However, the solution 
is not necessarily unique. If we take a solution (p and some a € that is constant on each connected 
component and vanishes on the boundary of T, it is clear that ^ + a is still an extension of ip and 


that 


dP{^) 

dip(ji) 


v=4> 


dP{^) 

dip{n) 


(p—cl)-\-(y 


so (p+a is another solution. We say that a connected component of a circuit touches the boundary 
if it contains a vertex in dN. Note that such an a must vanish on all connected components touching 
the boundary. 

With these preliminaries in hand, we can solve the Dirichlet problem: 


Proposition 2.6. For any boundary potential ip € there exists a potential cp obeying the 

principle of minimum power for ip. If we also demand that (p vanish on every connected component 
of r not touching the boundary, then (p is unique. 


Proof. For existence, observe that the power is a nonnegative quadratic form, the extensions of ip 
form an affine subspace of and a nonnegative quadratic form restricted to an affine subspace of 
a real vector space must reach a minimum somewhere on this subspace. 

For uniqueness, suppose that both (p and (p' obey the principle of minimum power for ip. Let 


a = (p' — (p. 


Then 


^ law ^\dN Ip 4’ 0, 

so 0 + Aa is an extension of ip for all A € M. This implies that 


/(A) := P{P + Xa) 


is a smooth function attaining its minimum value at both A = 0 and A = 1. In particular, this 
implies that f'{0) = 0. But this means that when writing / as a quadratic, the coefficient of A must 
be 0, so we can write 


2/(A) = + Aa)(t(e)) -{(p + Aa)(s(e)))^ 

= Y “ <^(s(e))) + A(a(i(e)) - a(s(e)))) 

= Y “ 4{s[e))f + A-term + A^ ^ -^(a(t(e)) - a(s(e)))^ 

= Y. " 4{.s{e))f + \^Y “ a(s(e)))"' 

eGE ^ ^ eGE ^ ' 
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Then 


/(I) - /(O) = -rT(«(^(e)) - «(s(e)))^ = 0, 

^ eGE 

so a{t{e)) = a(s(e)) for every edge e € E. This implies that a is constant on each connected 
component of the graph T of our circuit. 

Note that as a\dN = 0, a vanishes on every connected component of T touching the boundary. 
Note also that there always exists a solution (j) that vanishes on any connected component of T not 
touching the boundary: indeed, the total power dissipated is simply the sum of the power dissipated 
on each connected component, and any such vanishing potential dissipates zero power away from the 
boundary. Restricting our attention to solutions (j) and 4>' with this property, we see that a = (f>' — (f> 
vanishes on all connected components of T, and hence is identically zero. Thus (f>' = (j>, and this 
extra condition ensures a unique solution to the Dirichlet problem. □ 

We have also shown the following: 

Proposition 2.7. Suppose ip S and cp is a potential obeying the principle of minimum power 
for Ip. Then (p' obeys the principle of minimum power for ip if and only if the difference p' — p is 
constant on every connected component of T and vanishes on every connected component touching 
the boundary o/T. 

Furthermore, p depends linearly on p: 

Proposition 2.8. Fix p G , and suppose p S R'^ is the unique potential obeying the principle 
of minimum power for p that vanishes on all connected components of T not touching the boundary. 
Then p depends linearly on p. 

Proof. Fix p,p' G , and suppose p,p' G R^ obey the principle of minimum power for p,p' 
respectively, and that both p and p' vanish on all connected components of F not touching the 
boundary. 

Then, for all A G R, 


iP “b ) Ijjaw — + ^P Irsn — p A'0 

and 

(V((/l + A(/l'))|j^g„ = + A(V(/i')|Ra]v = 0. 

Thus p + Xp' solves the Dirichlet problem for p + Xp', and thus p depends linearly on p. □ 

Bamberg and Sternberg [8] describe another way to solve the Dirichlet problem, going back to 
Weyl [46]. 


2.4 Equivalent circuits 

We have seen that boundary potentials determine, essentially uniquely, the value of all the elec¬ 
tric properties across the entire circuit. But from the perspective of control theory, this internal 
structure is irrelevant: we can only access the circuit at its terminals, and hence only need con¬ 
cern ourselves with the relationship between boundary potentials and boundary currents. In this 
section we streamline our investigations above to state the precise way in which boundary currents 
depend on boundary potentials. In particular, we shall see that the relationship is completely cap¬ 
tured by the functional taking boundary potentials to the minimum power used by any extension of 
that boundary potential. Furthermore, each such power functional determines a different boundary 
potential-boundary current relationship, and so we can conclude that two circuits are equivalent if 
and only if they have the same power functional. 
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An ‘external behavior’, or behavior for short, is an equivalence class of circuits, where two are 
considered equivalent when the boundary current is the same function of the boundary potential. 
The idea is that the boundary current and boundary potential are all that can be observed ‘from 
outside’, i.e. by making measurements at the terminals. Restricting our attention to what can be 
observed by making measurements at the terminals amounts to treating a circuit as a ‘black box’: 
that is, treating its interior as hidden from view. So, two circuits give the same behavior when they 
behave the same as ‘black boxes’. 

First let us check that the boundary current is a function of the boundary potential. For this we 
introduce an important quadratic form on the space of boundary potentials: 

Definition 2.9. The power functional Q : R o/ a circuit with extended power functional 

P is given by 

Q(V’) = min P((/i). 

<l>\jtdN=V 

Proposition 2.6 shows the minimum above exists, so the power functional is well defined. Thanks 
to the principle of minimum power, Q{ip) equals ^ times the power dissipated by the circuit when 
the boundary voltage is ip. We will later see that in fact Q('!/') is a nonnegative quadratic form on 

Since Q is a smooth real-valued function on its differential dQ at any given point ijj G R^-^ 
defines an element of the dual space (R^'^)*, which we denote by dQ^. In fact, this element is equal 
to the boundary current b corresponding to the boundary voltage f): 

Proposition 2.10. Suppose tjj G R^'^. Suppose (j) is any extension of ijj minimizing the power. 
Then dQ^ G (R^^)* = R^-^ gives the boundary current of the current induced by the potential (f. 

Proof. Note first that while there may be several choices of f minimizing the power subject to the 
constraint that = ip, Proposition 2.7 says that the difference between any two choices vanishes 

on all components touching the boundary of F. Thus, these two choices give the same value for the 
boundary current i: dN —>■ R. So, with no loss of generality we may assume (p is the unique choice 
that vanishes on all components not touching the boundary. Write I: R for the extension of 

l: on to N taking value 0 on N \ ON. 

By Proposition 2.8, there is a linear operator 

/: ^ rN 

sending ip G R^^ to this choice of cp, and then 

Q{iP) = PifiP). 


16 


Given any ip' € 


we thus have 


dQ^{i;’) = —Q{cP + Xi;’) 


A=0 


= ^P{f{^ + X^P')) 


1 d 1 

2 dX r{e) 

eeB ^ ' 


A=0 


A=0 


fi'tp + X->p')){t{e)) - (/(V' + Xp}')){s{e))^ 

f'ip{s(e))) + A(/i/''(t(e)) -/V^'(s(e)))^ 


eeB 

1 


A=0 


= J 2 “ /V’(s(e)))(/V’'(^(e)) - fip'isie))) 

'r(e) 


e^E 


= X! -^(e)(/^'(i(e)) - /V''(s(e))) 

eg-E 

= XI ( XI X 

n^N Yi(e)=n s(e)—n J 

= X 

nGN 

= X 

n^dN 

This shows that dQ'p = i, as claimed. Note that this calculation explains why we inserted a factor 
of i in the definition of P: it cancels the factor of 2 obtained from differentiating a square. □ 

Note this only depends on Q, which makes no mention of the potentials at nonterminals. This is 
amazing: the way power depends on boundary potentials completely characterizes the way boundary 
currents depend on boundary potentials. In particular, in Part III we shall see that this allows us 
to define a composition rule for behaviors of circuits. 

To demonstrate these notions, we give a basic example of equivalent circuits. 

Example 2.11 (Resistors in series). Resistors are said to be placed in series if they are placed end to 
end or, more precisely, if they form a path with no self-intersections. It is well known that resistors 
in series are equivalent to a single resistor with resistance equal to the sum of their resistances. To 
prove this, consider the following circuit comprising two resistors in series, with input A and output 
C: 

fAB B ^BC 

A • -^^-• C 

Now, the extended power functional P: —>■ R for this circuit is 


PW = 




2 \rAB 

while the power functional Q: 


{cP{A) - p>{B)y 


tbc 


(cPiB) - cPiC))- 


is given by minimization over values of (j){B) = x: 


Qitp) = min 


(. 


1 


(tpiA) - x)^ 


1 


■{x-ipiC))' 


a:gR 2 \rAB ' ' TbC 

Differentiating with respect to a;, we see that this minimum occurs when 

— {x - ipiA)) + — {x - tp{C)) = 0, 
rAB Tbc 
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and hence when x is the r-weighted average of '0(A) and ■0(C'): 


TBC'ipiA) + TABi’iC) 

^ ^ 

rBC +tab 

Substituting this value for x into the expression for Q above and simplifying gives 

1 1 .. 


Q(0) = 


2 TAB +rBC 
This is also the power functional of the circuit 

rAB + rBC 


(0(A)-0(0)^ 


-• C 


and so the circuits are equivalent. 

2.5 Dirichlet forms 

In the previous subsection we claimed that power functionals are quadratic forms on the boundary 
of the circuit whose behavior they represent. They comprise, in fact, precisely those quadratic forms 
known as Dirichlet forms. 

Definition 2.12. Given a finite set S, a Dirichlet form on S is a quadratic form Q : R'® —>■ R 
given hy the formula 

<9(0) = 

i,j€S 

for some nonnegative real numbers Cij, and where we have written ipi = ip{i) G R. 

Note that we may assume without loss of generality that Cu = 0 and Cy = CjC, we do this hence¬ 
forth. Any Dirichlet form is nonnegative: Q(0) > 0 for all 0 G R'^. However, not all nonnegative 
quadratic forms are Dirichlet forms. For example, if S' = {1,2}, the nonnegative quadratic form 
Q(0) = (V'l + 02)^ is not a Dirichlet form. That said, the concept of Dirichlet form is vastly more 
general than the above definition: such quadratic forms are studied not just on finite-dimensional 
vector spaces R"® but on of any measure space. When this measure space is just a finite set, the 
concept of Dirichlet form reduces to the definition above. For a thorough introduction to Dirichlet 
forms, see the text by Fukushima [23]. For a fun tour of the underlying ideas, see the paper by 
Doyle and Snell [18]. 

The following characterizations of Dirichlet forms help illuminate the concept: 

Proposition 2.13. Given a finite set S and a quadratic form Q: R'® —>■ R, the following are 
equivalent: 

(i) Q is a Dirichlet form. 

(ii) Qif) < Q(0) whenever |0i — 0^1 < |0i — 0:;j for all i,j. 

(Hi) Q(0) = 0 whenever 0j is independent of i, and Q obeys the Markov property.' Q{<fi) < Q(0) 
when (f>i = min(0j, 1). 

Proof. See Fukushima [23]. □ 

While the extended power functionals of circuits are evidently Dirichlet forms, it is not immediate 
that all power functionals are. For this it is crucial that the property of being a Dirichlet form is 
preserved under minimising over linear subspaces of the domain that are generated by subsets of the 
given finite set. 


18 






Proposition 2.14. If Q: —)• K is a Dirichlet form, then 

min Q(—, v ): K'® ^ ffi. 


is Dirichlet. 


Proof. We first note that minj^gRS (5(—, z/) is a quadratic form. Again, minj^ggr (5(—, z/) is well 
defined as a nonnegative quadratic form also attains its minimum on an affine subspace of its 
domain. Furthermore min^^ggr Q(—, v) is itself a quadratic form, as the partial derivatives of Q are 
linear, and hence the points at which these minima are attained depend linearly on the argument of 
inin^^gRT Q{-,v). 

Now by Proposition 2.13, Q{(f) < Q{(f') whenever \4>i — (/zjl < |(/)- — </>'| for all i,j G S -\-T. In 
particular, this implies min^^gRT Q{ip,v) < min^^gRT Q{ip',v) whenever j-z/zi — t/’j l < It/’i — V'jl for 
i,j G S. Using Proposition 2.13 again then implies that min^^gRT Q{—, v) is a Dirichlet form. □ 


Corollary 2.15. Let Q : —^M.be the power functional for some circuit. Then Q is a Dirichlet 

form. 


Proof. The extended power functional P is a Dirichlet form, and writing © R^\^^ allows 

us to write 


Qi-) 


min 




□ 


The converse is also true: simply construct the circuit with set of vertices dN and an edge of 
resistance 2 ~ between any i,j€ dN such that the term Cijlppi — ipj) appears in the Dirichlet form. 
This gives: 

Proposition 2.16. A funetion Q : R^-^ —> R is the power funetional for some cireuit if and only if 
it is a Diriehlet form. 


This is an expression of the ‘star-mesh transform’, a well-known fact of electrical engineering 
stating that every circuit of linear resistors is equivalent to some complete graph of resistors between 
its terminals. For more details see [31]. We may interpret the proof of Proposition 2.14 as showing 
that intermediate potentials at minima depend linearly on boundary potentials, in fact a weighted 
average, and that substituting these into a quadratic form still gives a quadratic form. 


In summary, in this section we have shown the existence of a surjective function 


r circuits of linear resistors 
I with boundary dN 


Dirichlet forms on dN 


mapping two circuits to the same Dirichlet form if and only if they have the same external behavior. 
In the next section we extend this result to encompass inductors and capacitors too. 


3 Inductors and capacitors 

The intuition gleaned from the study of resistors carries over to inductors and capacitors too, to 
provide a framework for studying what are known as passive linear networks. To understand in¬ 
ductors and capacitors in this way, however, we must introduce a notion of time dependency and 
subsequently the Laplace transform, which allows us to work in the so-called frequency domain. 
Here, like resistors, inductors and capacitors simply impose a relationship of proportionality be¬ 
tween the voltages and currents that run across them. The constant of proportionality is known as 
the impedance of the component. 
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As for resistors, the interconnection of such components may be understood, at least formally, as 
a minimization of some quantity, and we may represent the behaviors of this class of circuits with a 
more general idea of Dirichlet form. We conclude this section by noting an obstruction to building 
a composition rule for Dirichlet forms, motivating our work in Part II. 


3.1 The frequency domain and Ohm’s law revisited 

In broadening the class of electrical circuit components under examination, we find ourselves dealing 
with components whose behaviors depend on the rates of change of current and voltage with respect 
to time. We thus now consider time-varying voltages v: [0,oo) —>■ K and currents i: [0,oo) —>■ R., 
where t G [0, oo) is a real variable representing time. For mathematical reasons, we restrict these 
voltages and currents to only those with (i) zero initial conditions (that is, /(O) = 0) and (ii) Laplace 
transform lying in the field 


R(s) = |z(s) 



P, Q polynomials over R in s. 



of real rational functions of one variable. While it is possible that physical voltages and currents 
might vary with time in a more general way, we restrict to these cases as the rational functions are, 
crucially, well behaved enough to form a field, and yet still general enough to provide arbitrarily 
close approximations to currents and voltages found in standard applications. 

An inductor is a two-terminal circuit component across which the voltage is proportional to 
the rate of change of the current. By convention we draw this as follows, with the inductance L the 
constant of proportionality:^ 

L 


Writing v^it) and ihit) for the voltage and current over time t across this component respectively, 
and using a dot to denote the derivative with respect to time t, we thus have the relationship 

VL(t) = LiLit). 

Permuting the roles of current and voltage, a capacitor is a two-terminal circuit component across 
which the current is proportional to the rate of change of the voltage. We draw this as follows, with 
the capacitance C the constant of proportionality: 

C 

*- 1 1 -* 

Writing vc{t), ic{t) for the voltage and current across the capacitor, this gives the equation 

ic{t) = C vc{t). 

We assume here that inductances L and capacitances C are positive real numbers. 

Although inductors and capacitors impose a linear relationship if we involve the derivatives of 
current and voltage, to mimic the above work on resistors we wish to have a constant of proportion¬ 
ality between functions representing the current and voltage themselves. Various integral transforms 
perform just this role; electrical engineers typically use the Laplace transform. This lets us write a 
function of time t instead as a function of frequencies s, and in doing so turns differentiation with 
respect to t into multiplication by s, and integration with respect to t into division by s. 

®We follow the standard convention of denoting inductance by the letter L, after the work of Heinrich Lenz and 
to avoid confusion with the I used for current. 
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In detail, given a function f{t): [0, oo) —)• R, we define the Laplace transform of / 

OO 

^{f}is) = j fit)e-^*dt. 

0 

We also use the notation £{/}(s) = F{s), denoting the Laplace transform of a function in upper 
case, and refer to the Laplace transforms as lying in the frequency domain or s-domain. For us, 
the three crucial properties of the Laplace transform are then: 

(i) linearity: £{a/ + bg}{s) = aF{s) + bG{s) for a, 6 e M; 

(ii) differentiation: £{/}(s) = sF{s) — /(O); 

(iii) integration: if g{t) = J* f{T)dT then G(s) = ^F{s). 

Writing F(s) and I{s) for the Laplace transform of the voltage v{t) and current i{t) across a compo¬ 
nent respectively, and recalling that by assumption v(t) = i(t) =0 for t < 0, the s-domain behaviors 
of components become, for a resistor of resistance R: 

F(s) = RI{s), 

for an inductor of inductance L: 

V (s) = sLI{s), 

and for a capacitor of capacitance C: 

Vis) = ^/(s). 

Note that for each component the voltage equals the current times a rational function of the real 
variable s, called the impedance and in general denoted by Z. Note also that the impedance is a 
positive real function, meaning that it lies in the set 

{Z € R(s) : Vs e C Re(s) > 0 ==^ Re(Z(s)) > 0}. 

While Z is a quotient of polynomials with real cofficients, in this definition we are applying it to 
complex values of s, and demanding that its real part be positive in the open left half-plane. Positive 
real functions were introduced by Otto Brune in 1931, and they play a basic role in circuit theory 
[ 12 ]. 

Indeed, Brune convincingly argued that for any conceivable passive linear component with two 
terminals we have this generalization of Ohm’s law: 

y(s) = z(s)i(s) 

where I S R(s) is the current, V € R.(s) is the voltage and the positive real function Z is the 
impedance of the component. As we shall see, generalizing from circuits of linear resistors to 
arbitrary passive linear circuits is just a matter of formally replacing resistances by impedances. 

As we consider passive linear circuits with more than two terminals, however, the coefhcients of 
our power functionals lie not just in the set of positive real functions, but in its closure ]R(s)+ under 
addition, multiplication, and division. Note that this closure is strictly smaller than M(s), as each 
positive real function is takes strictly positive real values on the positive real axis, so no matter how 
we take sums, products, and quotients of positive real functions it never results in the zero function 
0 G R(s). To generalise from linear resistors to passive linear circuits, we then replace the field R by 
the larger field R(s), and replace the set of positive reals, R+ = (0, oo), by the set R(s)+. From a 
mathematical perspective we might as well work with any field with a notion of ‘positive element’. 
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Definition 3.1. Given a field F, we define a set of positive elements for F to be any subset 
F'*' C F not containing 0 and closed under addition, multiplication, and division. 

Our first motivating example arises from circuits made of resistors. Here F = K is the field of 
real numbers and we take F+ = (0, oo). Our second motivating example arises from general passive 
linear circuits. Here F = ]R(s) is the field of rational functions in one real variable, and we take 
F+ = ]R(s)+ to be the closure of the positive real functions, as defined above. 

In all that follows, we fix a field F equipped with a set of positive elements F+. By a ‘circuit’, 
we shall henceforth mean a circuit over F’*', as explained in Definition 2.2. To fix the notation: 

Definition 3.2. A (passive linear) circuit is a graph s,t: E N with E as its set o/ edges 
and N as its set of nodes, equipped with function Z: E ^ F’*' assigning each edge an impedance, 
together with finite sets X, Y, and functions i: X ^ N and o: Y ^ N. We call the sets i{X), 
o(Y), and ON = i{X) U o{Y) the inputs, outputs, and terminals or boundary of the circuit, 
respectively. 

With these new components in hand, we are in a position to clarify the analogy between elec¬ 
tronics and mechanics. 

3.2 The mechanical analogy 

Now that we have introduced inductors and capacitors, it is worth taking another glance at the 
analogy chart in Section 1. What are the analogues of resistance, inductance and capacitance in 
mechanics? If we restrict attention to systems with translational degrees of freedom, the answer is 
given in the following chart. 


Electronics 

Mechanics (translation) 

charge Q 

position q 

current i = Q 

velocity v = q 

flux linkage A 

momentum p 

voltage V = X 

force F = p 

resistance R 

damping coefficient c 

inductance L 

mass m 

inverse capacitance C~^ 

spring constant k 


A famous example concerns an electric circuit with a resistor of resistance R, an inductor of 
inductance L, and a capacitor of capacitance C, all in series: 



We saw in Example 2.11 that for resistors in series, the resistances add. The same fact holds more 
generally for passive linear circuits, so the impedance of this circuit is the sum 

Z = sLY R+{sC)-^. 

Thus, the voltage across this circuit is related to the current through the circuit by 

H(s) = {sL + R+{sC)-^)I{s) 

If vff) and i{f) are the voltage and current as functions of time, we conclude that 


v{t) = L—i{t) + Ri{f) + C ^ 


i{s) ds 


0 
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It follows that 


LQ + RQ + C-^Q = V 

where Q{t) = J*i{t)ds has units of charge. As the chart above suggests, this equation is analogous 
to that of a damped harmonic oscillator: 


mq + cq + kq = F 


where m is the mass of the oscillator, c is the damping coefficient, k is the spring constant and F is 
a time-dependent external force. 

For details, and many more analogies of this sort, see the book by Karnopp, Margolis and 
Rosenberg [28] or Brown’s enormous text [11]. While it would be a distraction to discuss them 
further here, these analogies mean that our work applies to a wide class of networked systems, not 
just electrical circuits. 


3.3 Generalized Dirichlet forms 


To understand the behavior of passive linear circuits we need to understand how the behaviors 
of individual components, governed by Ohm’s law, fit together to give the behavior of an entire 
network. Kirchhoff’s laws still hold, and so does a version of the principle of minimum power. 

As before, to each passive linear circuit we associate a generalised Dirichlet form. Recall that we 
have fixed a field F equipped with a set F+ of positive elements (see Definition 3.1). 


Definition 3.3. Given a finite set S, a Dirichlet form over F on S' zs a quadratic form Q: F'® —)• F 
given hy the formula 

Qi'4’) ~ 'y ^ ^ijidl^i ~ '^o') 1 

ijes 

for some choice of Cij € F+ U {0}. 


Generalizing from circuits of resistors, we define the extended power functional P: F^ —>■ F 
of any circuit by 

Pi^p) = “ Pis{e))Y. 


2 ^ Z(e) 

e&E '• ' 


and we call ip G F^ a potential. Note that a held of characteristic 2 cannot be given a set of 
positive elements, so dividing by 2 is allowed. Note also that the extended power functional is a 
Dirichlet form on iV. 

Although it is not clear what it means to minimize over the held F, we can use formal derivatives 
to formulate an analogue of the principle of minimum power. This will actually be a ‘variational 
principle’, saying the derivative of the power functional vanishes with respect to certain variations in 
the potential. As before we shall see that given Ohm’s law, this principle is equivalent to Kirchhoff’s 
current law. 

Indeed, the extended power functional F((p) can be considered an element of the polynomial ring 
F[{(^(n)}neAr] generated by formal variables ip(n) corresponding to potentials at the nodes n G N. 
We may thus take formal derivatives of the extended power functional with respect to the p{n). 
We then call (f G F^ a realizable potential for the given circuit if for each nonterminal node, 
n G N \ dN, the formal partial derivative of the extended power functional with respect to (p(n} 
equals zero when evaluated at (p: 


dP 


d(p{n) 


= 0 




This terminology arises from the following fact, a generalization of Proposition 2.5: 
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Theorem 3.4. The potential (p G is a realizable potential for a given circuit if and only if the 
induced current 


= -yrriHtie)) - (/'(s(e))) 
Z(e) 


obeys Kirchhoff’s current law: 


Y, /(e) = E '<') 

s{e)—n t{e)—n 

for all n G N \ ON. 

Proof. The proof of this statement is exactly that for Proposition 2.5. 


□ 


A corollary of Theorem 3.4 is that the set of states—that is, potential-current pairs—that are 
compatible with the governing laws of a circuit is given by the set of realizable potentials together 
with their induced currents. 


3.4 A generalized minimizability result 

We begin to move from a discussion of the intrinsic behaviors of circuits to a discussion of their 
behaviors under composition. The key fact for composition of generalized Dirichlet forms is that, in 
analogy with Proposition 2.14, we may speak of a formal version of minimization of Dirichlet forms. 
We detail this here. In what follows let P be a Dirichlet form over F on some finite set S. 

Recall that given R G S, we call -0 € F'^ an extension of 0 S F^ if ip restricted to R equals ip. 
We call such an extension realizable if 


dP 

dip{s) 


= 0 

Lp = {jl 


for all s G 5* \ P. Note that over the real numbers R this means that among all the extensions of ip, 
Ip minimizes the function P. 


Theorem 3.5. Let P be a Dirichlet form over F on S, and let RG S be an inclusion of finite sets. 
Then we may uniquely define a Dirichlet form 

min P : F^ —)■ F 
S\R 


on R by sending each ip G F^ to the value P{ip) of any realizable extension ip of ip. 

To prove this theorem, we must first show that min^^^P is well defined as a function. 

Lemma 3.6. Let P be a Dirichlet form over F on S, let R G S be an inclusion of finite sets, and 
let Ip G F^. Then for all realizable extensions ip, ip' G F'® of ip we have P{ip) = P{ip'). 

Proof. This follows from the formal version of the multivariable Taylor theorem for polynomial 
rings over a field of characteristic zero. Let ip, ip' G F'^ be realizable extensions of ip, and note that 
dP^{ip — Ip') = 0, since for all s G P we have ip{s) — ip'{s) = 0, and for all s G P \ P we have 


dP 

d<p{s) 


= 0 . 


We may take the Taylor expansion of P around ip and evaluate at ip'. As P is a quadratic form, 
this gives 

P(0O = P(0) + dP^iip' - 0) + P(0' - 0) 

= P(0) + P{0' - 0). 
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Similarly, we arrive at 


P(V') = P(V'') + ^W-V'')- 

But again as P is a quadratic form, we then see that 

P{^') - P{^) = P{^' - -0) = P{^ - ^') = Pi-ip) - Pi-ip'). 

This implies that P(V'0 ~ -f’(V’) = Oj Eis required. 

It remains to show that min P remains a Dirichlet form. We do this inductively. 


□ 


Lemma 3.7. Let P be a Dirichlet form over F on S, and let s € S be an element of S. Then the 
map minj^j. P : F sending ip to Pi-ip) is a Dirichlet form on S \ {s}. 

Proof. Write j Cijiepi — (pj)"^- Observe that since i(pi — (pj)'^ = i(pj — (pi)'^, the definition 

of a Dirichlet form depends only on the sums Cy + Cji, and not intrinsically on the values of Cy 
and Cji. Previously, when defining the power functional of circuit, we chose coefficients such that 
Cy = Cy. Here, on the contrary and without loss of generality, it will be convenient for us to assume 
that Csk = 0 for all k. 

We then have 

dP ^ 

~ ^ (pk)-) 

v=4> k 


dps) 


and this is equal to zero when 




Ys — 

Xfc 

Note that the Cy lie in F+, and F+ is closed under addition, so X/c Cfcs 0. Thus min^gj P may be 
given explicitly by the expression 


min Pip = Cijip-p)‘^+ Y 

{s} 


*.ieS\{s} 


^GS\{s} 


Cis -ipi - 


Tdk ^ks-^k 
'P.k ^ka 


We must show this is a Dirichlet form on S' \ {s}. 

As the sum of Dirichlet forms is evidently Dirichlet, it suffices to check that the expression 


E 


Cis -ipi - 


Efc , 


is Dirichlet on S \ {s}. Multiplying through by the constant (Xfc Cks)'^ G F+, it further suffices to 
check 


Y^^‘> C-ks'4^i ^ ^ C-ks'4^k Cis T. CksiP - -tpk) 


I 


= E 




\ 


2 Yj ^ksCmsilpi - -ipk)iP - -tpm) + Y ^ksiP “ 


. fc,m 
y k^m 


) 




k,£,m 

k^m 


kl 
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is Dirichlet. But 


(V’fc - V'^)('0fe - '0m) + {fpe - - 0m) + (0m “ 0fe)(0m “ 0f) 

= + 0| + 0^ - i>k^i - 4’ki’m - 0^0m 

= 5((0fc - 0t)^ + (0fe - 0m)^ + (0<! - 0m)^), 


SO this expression is indeed Dirichlet. Indeed, pasting these computations together shows that 


niinP(0) = 
{*} 


E 


Cij + 


^is^js 

SfcCfes 


(0i - 0j)^- 


□ 


With these two lemmas, the proof of Theorem 3.5 becomes straightforward. 

Proof of Theorem 3.5. Lemma 3.6 shows that mingy/jP is a well-defined function. As P is a finite 
set, we may write it P = {si,..., s„} for some natural number n. Then we may define a sequence 
of functions Pi = min^^j Pi_i, 1 < i < n. Define also Pq = P, and note that P„ = min^^/j P. 
Then, by Lemma 3.7, each Pi is Dirichlet as Pi-i is. This proves the proposition. □ 

We can thus define the power functional of a circuit by analogy with circuits made of resistors: 

Definition 3.8. The power functional Q : R of a circuit with extended power functional 

P is given by 

Q = min P. 

N\dN 

As before, we define two circuits to be equivalent if they have the same power functional, and 
define the behavior of a circuit to be its equivalence class. 


3.5 Composition of Dirichlet forms 

It would be nice to have a category in which circuits are morphisms, and a category in which Dirichlet 
forms are morphisms, such that the map sending a circuit to its behavior is a functor. Here we 
present a naive attempt to constructed the category with Dirichlet forms as morphisms, using the 
principle of minimum power to compose these morphisms. Unfortunately the proposed category 
does not include identity morphisms. However, it points in the right direction, and underlines the 
importance of the cospan formalism we then turn to develop. 

We can define a composition rule for Dirichlet forms that reflects composition of circuits. Given 
finite sets S and T, let S + T denote their disjoint union. Let D{S, T) be the set of Dirichlet forms 
on 5" -I- T. There is a way to compose these Dirichlet forms 

o: P(P, U) X D(5, T) ^ D{S, U) 
defined as follows. Given P G D{T, U) and Q G D{S, T), let 

(PoQ)(a,7) = nnnQ(a,/3) -f P(/3,7), 

where a G P‘®,7 G P^. This operation has a clear interpretation in terms of electrical circuits: the 
power used by the entire circuit is just the sum of the power used by its parts. 

It is immediate from Theorem 3.5 that this composition rule is well defined: the composite of 
two Dirichlet forms is again a Dirichlet form. Moreover, this composition is associative. However, 
it fails to provide the structure of a category, as there is typically no Dirichlet form I5 G S) 
playing the role of the identity for this composition. For an indication of why this is so, let {•} be a 
set with one element, and suppose that some Dirichlet form /(/3,7) = fc(/3 — 7)^ G D({»}, {•}) acts 
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as an identity on the right for this composition. Then for all Q{a, /3) = c(a — /3)^ G £)({•}, {•}), we 
must have 


ca^ = Q(a, 0) 

= (/oQ)(a,0) 

= mm (3 (q;, /3) + /(/3, 0) 

= min k{a — /3)^ + c/3^ 
/SeF 

_ kc 2 


where we have noted that minimizes k{a — /3)^ + c/3^ with respect to (3. But for any choice of 

fc G F this equality only holds when c = 0, so no such Dirichlet form exists. Note, however, that for 
k >> c we have ca^ w '^3^0?■> so Dirichlet forms with large values of k —corresponding to resistors 
with resistance close to zero—act as ‘approximate identities’. 

In this way we might interpret the identities we wish to introduce into this category as the 
behaviors of idealized components with zero resistance: perfectly conductive wires. Unfortunately, 
the power functional of a purely conductive wire is undefined: the formula for it involves division 
by zero. In real life, coming close to this situation leads to the disaster that electricians call a ‘short 
circuit’: a huge amount of power dissipated for even a small voltage. This is why we have fuses and 
circuit breakers. 

Nonetheless, we have most of the structure required for a category. A ‘category without identity 
morphisms’ is called a semicategory, so we see 

Proposition 3.9. There is a semicategory where: 

• the objects are finite sets, 

• a morphism from T to S is a Dirichlet form Q G D{S, T). 

• composition of morphisms is given by 

{R o Q)( 7 , a) = nun (5(7,/3) + i?(/3, a). 

We would like to make this into a category. One easy way to do this is to formally adjoin 
identity morphisms; this trick works for any semicategory. However, we obtain a better category if 
we include more morphisms: more behaviors corresponding to circuits made of perfectly conductive 
wires. As the expression for the extended power functional includes the reciprocals of impedances, 
such circuits cannot be expressed within the framework we have developed thus far. Indeed, for 
these idealized circuits there is no function taking boundary potentials to boundary currents: the 
vanishing impedance would imply that any difference in potentials at the boundary induces ‘infinite’ 
currents. To deal with this issue, we generalize Dirichlet forms to Lagrangian relations. First, 
however, we develop a category theoretic framework, based around decorated cospans, to dehne the 
category of circuits itself and understand its basic properties. 


Part II 

Categories of Circuits 

In this part we move our focus from the semantics of circuit diagrams to the syntax, addressing 
the question “How do we interact with circuit diagrams?”. Informally, the answer to this is that 
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we interact with them by connecting them to each other, perhaps after moving them into the right 
form by rotating or reflecting them, or by crossing or bending some of the wires. To formalize 
this, we adopt a category theoretic viewpoint, defining various dagger compact categories with 
circuits and their behaviors as morphisms. We claim a formal analysis of this structure, especially 
of the composition or connection of circuits, has been overlooked in analysis of circuits thus far. 
This part culminates in the definition of two important categories, the category Circ of circuit 
diagrams, and the category LagrRel containing all behaviors of circuits. We also develop the technical 
material required to appreciate the structure of these categories, and that aids understanding of the 
relationship between the two, to be addressed in Part III. 

4 Decorated cospans 

We begin this part with a technical section describing a general technique for developing composition 
rules for structures on finite sets. As we have seen, whether represented by circuit diagrams or 
Dirichlet forms, circuits can be described as structures on a finite set of nodes. While this provides 
a good classification of the different types of circuits that exist, it does not allow for discussion of 
their composition. In this section, however, we describe a method for taking (I) a description of 
a structure that can be placed on finite sets together with (2) a description of how this structure 
interacts with functions between these sets, and producing a category which describes composition 
of structures. This category is built as a cospan category, with the apex of the cospan describing 
some structure, such as a circuit, and the feet of the cospan describing possible interfaces to this 
structure. 

4.1 Cospan categories 

Recall that a cospan from A to T in a category C is an object C in C with a pair of morphisms 
f: X ^ C, g: Y C: 



X Y. 

We call X and Y the feet of the cospan, and call C as its apex. When C has pushouts, cospans 
may be composed using the pushout from the common foot: given cospans X C Y and 

Y C 4^ Z, their composite cospan is X P Z where P, i: C ^ P, and i': C —> P 
form the top half of the pushout square: 


A 


P 



Y 


Z. 
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A map of cospans is a morphism h: C —>• C" in C between the apices of two cospans X C -(^ Y 
and X C Y with the same feet, such that 

C 



commutes. Given a category C with pushouts, we may thus define a category Cospan(C) with 
objects the objects of C and morphisms isomorphism classes of cospans. We will often abuse our 
terminology and refer to cospans themselves as morphisms in some cospan category Cospan(C); we 
of course refer instead to the isomorphism class of the said cospan. Note that there always exists an 
embedding 

C '—Cospan(C), 

with the embedding functor taking any object of C to its corresponding object in Cospan(C), and 
taking any morphism f: X ^ Y to the cospan X —^ Y Y. For this reason we often treat C as 

a subcategory of Cospan(C). 

Cospan categories also come equipped with a so-called ‘dagger structure’, which maps a cospan 

X C Y to its reflection Y —^ C G— X. Moreover, when C has finite colimits, the category 
Cospan(C) is compact. Since these facts are important in circuit theory, they warrant a brief review 
of dagger compact categories. 

4.2 Dagger compact categories 

In a dagger category, the distinction between the source and target of a morphism is arbitrary, and 
can be reversed. This is true of electrical circuits: if we like we may treat the set of inputs as the 
set of outputs instead, and the set of outputs as the set of inputs. 

Recall that a dagger category is a category equipped with an involutive, contravariant endo- 
functor that is the identity on objects [1, 40]. In other words, a dagger structure on a category C 
is a contravariant functor (—)!: C ^ C such that = A for all objects A G ObC and (/^)^ = / for 
all morphisms / in C. A dagger category is then a category equipped with a dagger structure. 

When other structure is present, we prefer this dagger to be compatible with it. We say that a 
morphism / is unitary if its dagger provides it with an inverse morphism. A symmetric monoidal 
dagger category is a symmetric monoidal category equipped with a symmetric monoidal dag¬ 
ger structure —that is, a dagger structure that coherently preserves the symmetric monoidal struc¬ 
ture. Concretely, this requires that the functor (—: C —> C°p be symmetric monoidal, and that the 
associator, unitors, and braiding of C be unitary. Furthermore, letting L and R be dual objects of 
a symmetric monoidal dagger category, with monoidal unit I, braiding aL,R ■ L ® R ^ R® and 
unit rj: I ^ R® L and counit e: L ® R ^ I, we say that L and R are dagger dual if ry = cr o eb 
A dagger compact category is a symmetric monoidal dagger category in which every object has 
a dagger dual, while a dagger functor is a functor F between dagger compact categories that 
preserves the dagger structures: F((—)1) = (F(—))b 

Importantly for our applications, dagger compact categories come with a graphical calculus, 
where each morphism is represented by a diagram such that two diagrams are considered equal by 
the rules of this calculus if and only if they are equal according to the defining laws of dagger compact 
categories [40]. In brief, to set up our conventions, we represent a morphism /: Ai (8> • • • 0 Xn —>■ 
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Yi (g) • • • (g) Ym as a ‘downwards flow chart’: 







/ 





Vl Yr. 


Composition is then represented by connecting the lines (or wires) representing the codomain of one 
morphism with the domain of the another placed below it, the monoidal product of two morphisms 
is represented by their side-by-side juxtaposition, the swap map by crossing two wires, the compact 
structure by bending a wire 180 degrees, and the dagger functor by flipping a diagram in the 
horizontal axis. We believe these operations on diagrams—placing diagrams on the same page, 
rearranging their wires/flipping them, and then connecting their wires to form a larger diagram— 
represents the collection of operations used for reasoning with circuit diagrams, and hence that 
dagger compact categories are an appropriate structure for the formalization of such. 

Example 4.1. As we have already remarked upon, whenever C is a category with finite colimits, 
Cospan(C) becomes a dagger compact category. We explain this for the main example of interest, 
the category Cospan(FinSet) of cospans in FinSet, the category of finite sets and functions. 

The category Cospan(FinSet) has a dagger structure that maps any cospan X C Y to 

the cospan Y —^ C X. It becomes symmetric monoidal with the tensor product arising from 
coproducts in the category of finite sets and the symmetric monoidal structure maps inherited from 
viewing (FinSet, -I-) as a subcategory. It is compact with each finite set X dual to any isomorphic 
set X with unit and counit 





0 x-i-a: X + X 0 


respectively, where ! is the unique function of the given type, a;: X AT is an isomorphism, and we 
write [/, 5 ] for the coproduct of morphisms / and g. It is a simple computation to check that dagger 
structure interacts with the symmetric monoidal structure and the duals for objects to indeed define 
a dagger compact category. 

We will spend much time discussing categories built from Cospan(FinSet). In the next subsection 
we describe methods to construct such categories. 

4.3 Decorated cospan categories 

While cospans in FinSet provide a good language to describe connections between finite sets—we 
view the feet as indexing connection points of the apex—to provide a framework for describing 
composition of more interesting structures than finite sets we need a way of incorporating this extra 
data. This is provided by the idea of an ‘F-decorated’ cospan: a cospan in FinSet in which the 
apex N is equipped with an element 1 —>■ FX of some object FN. We think of F as describing the 
collection of available structures on X, with examples such as the collection of circuit diagrams or 
Dirichlet forms on X in mind. 

When F is a lax symmetric monoidal functor, such decorated cospans are the morphisms of a 
dagger compact category. We assure the reader that all constructions in this section are well-defined; 
the interested reader can find details in Fong [20]. 
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Definition 4.2. Let 


{F, ip) : (FinSet, +) —^ [V, 

be a lax monoidal functor. We may define a category i^Cospan, the category of i^-decorated 
cospans, with objects being finite sets and morphisms being equivalence classes of pairs 

{X ^ N ^Y, 1 ^ FN) 

comprising a cospan X ^ N Y in FinSet together with an element of FN ; that is, a morphism 
I ^ FN where I is the unit for the tensor product in F. We call 1 FN the decoration of the 
decorated cospan. The equivalence relation arises from isomorphism of cospans; an isomorphism of 
cospans induces a one-to-one correspondence between their decorations. 

Composition in this category is given via pushout of cospans in FinSet; 



paired with the pushforward 

^ FN ® FM F{N + M) f{N Yy M) 


of the tensor product of the decorations along the coproduct of the pushout maps. 

While at first glance it might seem surprising that we can construct a composition rule for 
decorations s: 1 ^ FN and t: 1 ^ FM just from monoidal structure, the copair [JnCm] '. N + M ^ 
N +Y M of the pushout maps contains the data necessary to compose them. Indeed, this is the 
key insight of the decorated cospans construction. To wit, the coherence maps for the lax monoidal 
functor allow us to construct an element of F{N + M) from the monoidal product s ® t oi the 
decorations, and we may then post-compose with Jm] to arrive at an element of F{N +y M). 
The map [iw, Jm] encodes the identification of the image of F in with the image of the same in 
M, and so describes the ‘connection’ of the two ‘circuits’. 

Decorated cospan categories are so named as they generalize the category of cospans of finite 
sets: there is a faithful and essentially surjective functor 


Cospan(FinSet) '—5> FCospan. 


Details may be found in [20]. 

When the functor F is lax symmetric monoidal, the category FCospan becomes a dagger compact 
category. To be more precise, we define the tensor product of objects X and Y of FCospan to be 
their disjoint union X + Y, and define the tensor product of decorated cospans {X N 
Y, 1 ^ FN) and {X' "-^N' t^Y', 1-4 FN') to be 

F{N + N') \ 

/ 

We also write + for the tensor product in FCospan. 


/ N yN' 



\ XYX' Y YY' 
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The dagger structure for FCospan mimics that for Cospan(FinSet); it reflects the cospan part 
of a decorated cospan, keeping the same decoration: 

]{X ^ N 1 ^ FN) = {Y ^ N ^ X, 1 ^ FN). 

This gives the following fact: 

Lemma 4.3. Let F be a lax symmetric monoidal functor. Then with the above structure, the 
category T’Cospan is a dagger compact category. 

Proof. This is a special case of Theorem 3.4 of [20]. □ 

4.4 Functors between decorated cospan categories 

Decorated cospans allow us to understand the diagrammatic nature of structures on hnite sets, 
such as circuit diagrams. Equally crucial to our understanding of circuit diagrams, however, is 
their semantics as discussed in the previous section, say via their interpretation as Dirichlet forms. 
Decorated cospan categories also give us the tools to understand these interpretations through dagger 
compact categories. A consequence of the following general result is that the process of interpreting 
of a circuit diagram as a Dirichlet form preserves this dagger compact structure. 

Lemma 4.4. Let 

(F t), (G, 7 ): (FinSet, +) —^ (C, 0 ) 
be lax symmetric monoidal functors and let 

9: {F,^)^{G,j) 

be a monoidal natural transformation between them. Then we may define a functor, in fact a strict 
symmetric monoidal dagger functor. 


T: FCospan —GCospan 

by letting any finite set X in ECospan map to the same finite set as an object of GCospan, and 
letting any morphism 

(A ^ TV y, 1 ^ FN) 

map to: 

{X ^ N ^Y,l GN). 

Proof. This is a special case of Theorem 4.1 of [20]. □ 


5 The category of passive linear circuits 


In Part I we defined a circuit of linear resistors to be a labelled graph with marked input and output 
terminals, as in the example: 
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We then defined general passive linear circuits by replacing resistances with impedances chosen from 
a set of positive elements F+ in any field F. In fact, these circuits are examples of decorated cospans. 
This gives a dagger compact category Circ whose morphisms are circuits, with the dagger compact 
structure expressing standard operations on circuits. 

We actually give two constructions of this category, first arriving at a category of cospans deco¬ 
rated by F+-graphs, and then showing that this is a full subcategory of the decategorification of a 
bicategory of cospans of F+-graphs. 

5.1 A decorated cospan construction 

We officially defined a ‘circuit’ in Definition 3.2, but now we can give an equivalent definition using 
cospans: 

Lemma 5.1. A circuit is a cospan of finite sets X —^ N Y together with an ¥'^-graph whose 
set of nodes is N. 

Proof. This is just a matter of remembering the terminology: recall that an F^-graph is a graph 
whose edges are labelled with elements of F+, our chosen set of positive elements in the field F. □ 

This suggests that circuits should be morphisms in a decorated cospan category. Indeed, we can 
show that the map taking a finite set N to the set of F+-graphs with set N of nodes in fact forms 
a lax symmetric monoidal functor. This allows us to apply Lemma 4.2 to construct a category of 
circuits. 

To this end, define the functor 

Circuit: (FinSet,-!-) —> (Set, x) 

to take a finite set TV, as an object of FinSet, to the set Circuit(TV) of F+-graphs (TV, E, s, t, r) with 
TV as their set of nodes. On morphisms let it take a function /: TV —> TIL to the function that pushes 
labelled graph structures on a set TV forward onto the set M: 

Circuit(/): Circuit(TV) —> Circuit(T\L); 

(TV, E, s, t, r) I— > {M,E,fos,fo t, r). 

Note that as this map simply acts by post-composition, our map Circuit is indeed functorial. 

We then arrive at a lax symmetric monoidal functor by equipping this functor with the natural 
transformation 

Pn,m' Circuit(TV) x Circuit(TlL) —Circuit(TV-|-TIL); 

((TV, E, s, t, r), {M, E, s', t', r')) \—^ [N -\- M, E -\- F, s -\- s',t -\-1', [r, r ']), 

together with the unit map 


Pi : 1 —Circuit(0); 

• ^ ( 0 , 0 , 0 , 0 , 0 ), 

where we use 0 to denote both the empty set and the unique function of the appropriate codomain 
with domain the empty set. The naturality of this collection of morphisms, as well as the coherence 
laws for lax symmetric monoidal functors, follow from the universal property of the coproduct. 

Definition 5.2. We define 

Circ = CircuitCospan. 

Lemma 4.3 then specialises to: 
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Corollary 5.3. The category Circ is a dagger compact category. 

The different structures of this category capture different operations that can be performed with 
circuits. The composition expresses the fact that we can connect the outputs of one circuit to the 
inputs of the next, while the monoidal composition models the placement of circuits side-by-side. 
The symmetric monoidal structure allows us reorder input and output wires, and the compactness 
captures the interchangeability between input and outputs of circuits—that is, the fact that we can 
choose any input to our circuit and consider it instead as an output, and vice versa. Finally, the 
dagger structure expresses the fact that we may reflect a whole circuit, switching all inputs with all 
outputs. 


5.2 A bicategory of circuits 

As shown by Benabou [9], cospans are most naturally thought of as 1-morphisms in a bicategory. 
To obtain the category we are calling Cospan(C), we ‘decategorify’ this bicategory by discarding 
2-morphisms and identifying isomorphic 1-morphisms. Since we are studying circuits using deco¬ 
rated cospans, this suggests that circuits, too, are most naturally thought of as 1-morphisms in a 
bicategory. Indeed this is the case. 

One route to the bicategory whose 1-morphisms are circuits would be to take the theory of 
decorated cospans [20] and show that it can be enhanced to create not merely categories whose 
morphisms are isomorphism classes of decorated cospans, but bicategories whose 1-morphisms are 
exactly decorated cospans. A less powerful but easier approach is as follows. 

Recall that we define a graph to be a pair of functions s,t: E ^ N where E and N are finite 
sets, and an L-graph to be a graph further equipped with a function r: E ^ L. Thus, an L-graph 
looks like this: 

L E N 

t 

Given L-graphs F = {E, N, s,t,r) and F' = {E',N',s',t',r'), a morphism of L-graphs F —>■ F' is a 
pair of functions e: E ^ E', v: N ^ N' such that the following diagrams commute: 



These L-graphs and their morphisms form a category L-Graph. Using results about colimits in the 
category of sets, it is straightforward to check that this category has finite colimits. 

Benabou [9] gave a general construction of bicategories starting from any category with pushouts. 
Since L-Graph has pushouts, it follows that there is a bicategory Cospan(L-Graph) with 

• L-labelled graphs as objects, 

• cospans in L-Graphas morphisms, and 

• maps of cospans as 2-morphisms. 

Since L-Graph also has coproducts, this bicategory is symmetric monoidal and in fact compact, 
thanks to the work of Stay [44]. We make the following definition: 

Definition 5.4. The bicategory 2-Circ is the full and 2-full sub-bicategory o/Gospan(F+-Graph) 
with objects those ¥'^ -graphs with no edges. 
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It can be shown that every object in Cospan(F+-Graph) is self-dual, and this implies that 2-Circ 
is again a compact closed bicategory. Note that the category obtained by decategorifying this 
bicategory has finite sets as objects, with morphisms being isomorphism classes of cospans in F+- 
Graph with feet such objects. Thus, decategorifying 2-Girc gives a category equivalent to our 
previously defined category Girc. 


6 Circuits as Lagrangian relations 

In the first part of this paper, we explored the semantic content contained in circuit diagrams, leading 
to an understanding of circuit diagrams as expressing some relationship between the potentials and 
currents that can simultaneously be imposed on some subset, the so-called terminals, of the nodes 
of the circuit. We called this collection of possible relationships the behavior of the circuit. While in 
that setting we used the concept of Dirichlet forms to describe this relationship, we saw in the end 
that describing circuits as Dirichlet forms does not allow for a straightforward notion of composition 
of circuits. 

In this section, inspired by the principle of least action of classical mechanics in analogy with the 
principle of minimum power, we develop a setting for describing behaviors that allows for easy discus¬ 
sion of composite behaviors: Lagrangian subspaces of symplectic vector spaces. These Lagrangian 
subspaces provide a more direct, invariant perspective, comprising precisely the set of vectors de¬ 
scribing the possible simultaneous potential and current readings at all terminals of a given circuit. 
As we shall see, one immediate and important advantage of this setting is that we may model wires 
of zero resistance. 

Recall that we write F for some field, which for our applications is usually the field R of real 
numbers or the field ]R(s) of rational functions of one real variable. 

6.1 Symplectic vector spaces 

A circuit made up of wires of positive resistance defines a function from boundary potentials to 
boundary currents. A wire of zero resistance, however, does not define a function: the principle of 
minimum power is obeyed as long as the potentials at the two ends of the wire are equal. More 
generally, we may thus think of circuits as specifying a set of allowed voltage-current pairs, or as 
a relation between boundary potentials and boundary currents. This set forms what is called a 
Lagrangian subspace, and is given by the graph of the differential of the power functional. More 
generally, Lagrangian submanifolds graph derivatives of smooth functions: they describe the point 
evaluated and the tangent to that point within the same space. 

The material in this section is all known, and follows without great difficulty from the definitions. 
To keep this section brief we omit proofs. See any introduction to symplectic vector spaces, such as 
Gimasoni and Turaev [16] or Piccione and Tausk [36], for details. 

Definition 6.1. Given a finite-dimensional vector space V over a field F, a symplectic form 
w: P X R — > F on V is an alternating nondegenerate bilinear form. That is, a symplectic form lo is 
a function P x R —^ F that is 

(i) bilinear: for all X g¥ and all u,v gV we have a;(Au,u) = uj{u,Xv) = Xuj{u,v); 

(ii) alternating: for all v G V we have u}(v,v) = 0; and 

(Hi) nondegenerate: given v gV, uj(u,v) = 0 for all u G V if and only ifu = 0. 

A symplectic vector space {V,uj) is a vector space V equipped with a symplectic form uj. 

Given symplectic vector spaces (Vi,a;i), (V 2 , W 2 ), a symplectic map is a linear map 

f: (Vi,wi) —(V2,W2) 
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such that ui 2 {f (u), f {v)) = uji{u,v) for all u,v €Vi. A symplectomorphism is a symplectic map 
that is also an isomorphism. 

An alternating form is always antisymmetric, meaning that uj{u, v) = —uj{v, u) for all u,v G V. 
The converse is true except in characteristic 2. A symplectic basis for a symplectic vector space 
(y, to) is a basis {pi,... ,pn, < 71 ,, qn} such that uj{pi,pj) = uj{qi, qj) = 0 for all 1 < *, j < n, and 
u!{pi,qj) = dij for all 1 < < n, where Sij is the Kronecker delta, equal to 1 when i = j, and 0 

otherwise. A symplectomorphism maps symplectic bases to symplectic bases, and conversely, any 
map that takes a symplectic basis to another symplectic basis is a symplectomorphism. 

Example 6.2 (The symplectic vector space generated by a finite set). Given a finite set N, we 
consider the vector space © (F-^) a symplectic vector space (F^ © (F^) ,w), with symplectic 
form 

w((^, f), i')) = i'{(l>) - 

Let {(j)n}n&N be the basis of F-^ consisting of the functions N ^ ¥ mapping n to 1 and all other 
elements of N to 0, and let {in}nGN C (F-^)* be the dual basis. Then {(^„,0), (0,z„)}„gAr forms a 
symplectic basis for F'^ © (F-^)*. 

There are two common ways we will build symplectic spaces from other symplectic spaces: con¬ 
jugation and summation. Given a symplectic form to on V, we may define its conjugate symplectic 
form to = —to, and write the conjugate symplectic space (y,Zo) as V. Given two symplectic vector 
spaces {U, v), {V, w), we consider their direct sum [/© V a symplectic vector space with the symplec¬ 
tic form V + OJ, and call this the sum of the two symplectic vector spaces. Note that this is neither 
the product nor coproduct in the category of symplectic vector spaces and symplectic maps. 

The symplectic form provides a notion of orthogonal complement. Given a subspace S of V , we 
define its complement 

S'° = {u G y I w(n, s) = 0 for all s G S}. 

Note that this construction obeys the following identities, where S and T are subspaces of V : 

dim S + dim S° = dim V 
(S°)° = S 
(S + T)° = s° nT° 

(SnT)° = S°+T°. 

In the symplectic vector space F^ © (F^)*, the subspace F^ has the property of being a maximal 
subspace such that the symplectic form restricts to the zero form on this subspace. Subspaces with 
this property are known as Lagrangian subspaces, and they may all be realized as the image of 
F'^ © (F^)* under symplectomorphisms from F^ © (F^)* to itself. 

Definition 6.3. Let S be a linear subspace of a symplectic vector space (y,uj). We say that S is 
isotropic ifuj\sxS = 0, and that S is coisotropic if S° is isotropic. A subspace is Lagrangian if 
it is both isotropic and coisotropic, or equivalently, if it is a maximal isotropic subspace. 

Lagrangian subspaces are also known as Lagrangian correspondences and canonical relations. 
Note that a subspace S is isotropic if and only if S' C S° . This fact helps with the following 
characterizations of Lagrangian subspaces. 

Proposition 6.4. Given a subspace L C V of a symplectic vector space (y, w), the following are 
equivalent: 

(i) L is Lagrangian. 

(a) L is maximally isotropic. 
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(in) L is minimally coisotropic. 

(iv) L = L°. 

(v) L is isotropic and dimL = ^dimV. 

From this proposition it follows easily that the direct sum of two Lagrangian subspaces is La- 
grangian in the sum of their ambient spaces. We also observe that an advantage of isotropy is that 
there is a good way to take a quotient of a symplectic vector space by an isotropic subspace—that 
is, there is a way to put a natural symplectic structure on the quotient space. 

Proposition 6.5. Let S be an isotropic subspace of a sympleetic vector space {V,u}). Then S°/S is 
a symplectic vector space with symplectic form uj'{v + 5”, u + 5) = ui{v, u). 

Proof. The function lo' is well defined due to the isotropy of S: by definition, adding any pair (s, s') 
of elements of S' to a pair {v, u) of elements of S° does not change the value of w(u + s, u + s'). As 
w is a symplectic form, one can check that uj' is too. □ 


6.2 Lagrangian subspaces from quadratic forms 

Lagrangian subspaces are of relevance to us here as the behavior of any passive linear circuit forms a 
Lagrangian subspace of the symplectic vector space generated by the nodes of the circuit. We think 
of this vector space as comprising two parts: a space of potentials at each node, and a dual space 
(F^) of currents. To make clear how circuits can be interpreted as Lagrangian subspaces, here 
we describe how Dirichlet forms on a finite set N give rise to Lagrangian subspaces of F^ © (F^) . 
More generally, we show that there is a one-to-one correspondence between Lagrangian subspaces 
and quadratic forms. 


Proposition 6.6. Let N be a finite set. Given a quadratie form Q over F on N, the subspace 

Lq = {icf,dQ^) I 0 G F"^} C F^ © (F^)*, 


where dQ^ G (F^)* is the formal differential of Q at (j) € F^, is Lagrangian. Moreover, this 
construction gives a one-to-one correspondence 


Quadratic forms over ¥ on N 


Lagrangian subspaces of¥^ © (F^)” 
with trivial intersection with 
{0}©(F^)* CF^©(F^)* 


Proof. The symplectic structure on F^ © (F^)* and our notation for it is given in Example 6.2. 
Note that for all n,m € N the corresponding basis elements 


d^Q 

bd((>rin 




SO dQ^{if) = dQ.tj,{(j)) for all </>, V’ € F^. Thus Lq is indeed Lagrangian: for all V' S F^ 

dQ,j)), (gp, dQ,p)') — dQ^pi^cp) dQp(ip) — 0 . 


Observe also that for all quadratic forms Q we have dQo = 0, so the only element of Lq of the 
form (0, i), where i G (F^)*, is (0, 0). Thus Lq has trivial intersection with the subspace {0}©(F^)* 
of F-^ © (F'^)*. This Lq construction forms the leftward direction of the above correspondence. 

For the rightward direction, suppose that L is a Lagrangian subspace of F'^ © (F'^)* such that 
L n ({0} © (F-^)*) = {(0,0)}. Then for each (p G F^^, there exists a unique ip G (F-^)* such that 
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(</>, € L. Indeed, if and were distinct elements of (F'^)* with this property, then by linearity 

{Ojifj, — i'^) would be a nonzero element of L fl ({0} © (F^)*), contradicting the hypothesis about 
trivial intersection. We thus can define a function, indeed a linear map, F^ —>■ (F-^)*; (p i—>■ This 
defines a bilinear form Q{(j),'tp) = on F^ ©F^, and so Q{4>) = defines a quadratic form 

on F^. 

Moreover, L is Lagrangian, so 

w(((/), i^), (V^, i'lp)) = = 0, 

and so Q{—,—) is a symmetric bilinear form. This gives a one-to-one correspondence between 
Lagrangian subspaces of specified type, symmetric bilinear forms, and quadratic forms, and so in 
particular gives the claimed one-to-one correspondence. □ 

In particular, every Dirichlet form defines a Lagrangian subspace. 

6.3 Lagrangian relations 

Recall that a relation between sets X and T is a subset R of their product X x Y. Furthermore, 
given relations R C X x Y and S C Y x Z, there is a composite relation {S o R) C X x Z given 
by pairs (x, z) such that there exists y € Y with (x, y) € R and {y, z) G S —a direct generalization 
of function composition. A Lagrangian relation between symplectic vector spaces Vi and V 2 is a 
relation between Vi and V 2 that forms a Lagrangian subspace of the symplectic vector space Vi © V 2 ■ 
This gives us a way to think of certain Lagrangian subspaces, such as those arising from circuits, as 
morphisms, giving a way to compose them. 

Definition 6.7. A Lagrangian relation L: Vi —>■ V 2 is a Lagrangian subspace L o/Vi © V 2 . 

This is a generalization of the notion of symplectomorphism: any symplectomorphism /: Vi —> V 2 
forms a Lagrangian subspace when viewed as a relation / C Vi © V 2 - More generally, any symplectic 
map /: Vi —V 2 forms an isotropic subspace when viewed as a relation in Vi © V 2 ■ 

Importantly for us, the composite of two Lagrangian relations is again a Lagrangian relation. 
This is well known [45], but sufficiently easy and important to us that we provide a proof. 

Proposition 6 . 8 . Let L: Vi —>■ V 2 and L': V 2 —be Lagrangian relations. Then their composite 
relation L' o L is a Lagrangian relation Vi —>■ V 3 . 

We prove this proposition by way of two lemmas detailing how the Lagrangian property is 
preserved under various operations. The first lemma says that the intersection of a Lagrangian 
space with a coisotropic space is in some sense Lagrangian, once we account for the complement. 

Lemma 6.9. Let L QV be a Lagrangian subspace of a symplectic vector space V, and S Q V be an 
isotropic subspace ofV. Then [L n S°) + S' C 1/ is Lagrangian in V. 

Proof. Recall from Proposition 6.4 that a subspace is Lagrangian if and only if it is equal to its 
complement. The lemma is then immediate from the way taking the symplectic complement interacts 
with sums and intersections: 

{{Lns°)+s)° = {Lns°)°ns° = (L°+(s°)°)ns° = (L+s)ns° = (Lns°)+(sns°) = (Lns°)+s. 

Since (L n S°) + S is equal to its complement, it is Lagrangian. □ 

The second lemma says that if a subspace of a coisotropic space is Lagrangian, taking quotients 
by the complementary isotropic space does not affect this. 
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Lemma 6.10. Let L QV he a Lagrangian subspaee of a symplectic vector space V, and S Q L an 
isotropic subspace ofV contained in L. Then L/S C S°/S is Lagrangian in the quotient symplectic 
space S°/S. 

Proof. As L is isotropic and the symplectic form on S°/S is given by uj'{v + S,u + S) = uj{v, u), the 
quotient L/S is immediately isotropic. Recall from Proposition 6.4 that an isotropic subspace S of 
a symplectic vector space V is Lagrangian if and only if dim -S' = ^ dim V. Also recall that for any 
subspace dim 5' + dim 5'° = dimP. Thus 


dim(L/S') = dim L — dim -S' = | dim V — dim S 

= ^(dimS" + dimS"®) — dim^” = i(dim5'° — dimS') = ^ dim(S'°/S'). 

Thus L/S is Lagrangian in S°/S. □ 

Combining these two lemmas gives a proof that the composite of two Lagrangian relations is 
again a Lagrangian relation. 

Proof of Proposition 6.8. Let A be the diagonal subspace 


A = {( 0 , V2, W2, 0 ) \v2 G V2} C Vi 0 V2 © V2 © V3. 

Observe that A is isotropic, and has coisotropic complement 

A° = {{vi,V2,V2,V^) I Vi G Vi} c 14 © V 2 © V 2 © V 3 . 

As A is the kernel of the restriction of the projection map Vi © V 2 © V 2 © VJs — >■ Fi © 14 to A°, and 
after restriction this map is still surjective, the quotient space A°/A is isomorphic to Vi © V 3 . 

Now, by definition of composition of relations, 

L' o L = {(ni,r! 3 ) | there exists V2 G V 2 such that {vi,V2) G L, {v2,V3) G L'}. 


But note also that 

L®L' = {{vi,V2,V2,V3) I (■!;i,'!;2 ) G L, {v2,V3) G L'}, 


so 


(L © L') n A° = {{vi,V2,V2,V3) I there exists V2 G V2 such that (^1,^2) G L, {v2,V3) G L'}. 
Quotienting by A then gives 


L' o L = ((L © L') n A°) + A)/A. 

As L'(BL is Lagrangian in Vi©V 2 ©^©V 3 , Lemma 6.9 says that (L'©L)nA°) + A is also Lagrangian 
in Vi © V 2 © V 2 © V 3 . Lemma 6.10 thus shows that L' o L is Lagrangian in A°/A = 14 © 14, as 
required. □ 

Note that this composition is associative. We shall see later that this composition agrees with 
composition of Dirichlet forms, and hence also composition of circuits. 
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6.4 The dagger compact category of Lagrangian relations 

Lagrangian relations solve the identity problems we had with Dirichlet forms: given a symplectic 
vector space V, the Lagrangian relation id: V ^ V specified by the Lagrangian subspace 

id={{v,v) \ v €V} CV ®V, 

acts as an identity for composition of relations. We thus have a category. 

Definition 6.11. We write LagrRel for the category with symplectic vector spaces as objects and 
Lagrangian relations as morphisms. 

In fact the move to the setting of Lagrangian relations, rather than Dirichlet forms, adds far 
richer structure than just identity morphisms. The category LagrRel can be viewed as endowed 
with the structure of a dagger compact category. We lay this out in steps. 

Symmetric monoidal structure 

We define the tensor product of two objects of LagrRel to be their direct sum. Similarly, we define 
the tensor product of two morphisms L:U^V,LCU(BV and K: T ^ W, RT C T © W to be 
their direct sum 

L®K CU®V®T®W, 

but considered as a subspace of the naturally isomorphic space U (BT(BV(BW. Despite this subtlety, 
we abuse our notation and write their tensor product L® K: U ®T —^-R© W, and move on having 
sounded this note of caution. 

Note that the direct sum of two Lagrangian subspaces is again Lagrangian in the direct sum of 
their ambient spaces, and the zero dimensional vector space { 0 } acts as an identity for direct sum. 
Indeed, defining for all objects U, V, W in LagrRel unitors: 

Av = {(0,u,u)}C{0}©R©R, 

Pv = {(v, 0, u)} C R © {0} © R, 


associators: 

ctu,v,w = {(w, V, w, u, V, lu)} C ([/ © R) © IR © [/ © (R © W), 

and braidings: 

= {(u, v,v,u) \ u & U,v ffU ®V ®V ®U, 

we have a symmetric monoidal category. Note that all these structure maps come from symplec- 
tomorphisms between the domain and codomain. From this viewpoint it is immediate that all the 
necessary diagrams commute, so we have a symmetric monoidal category. 

Duals for objects 

Each object R of LagrRel is dual to its conjugate space R, with cup ?7:{0 }—>'R©R given by 

?7 = {(0, V, v) I V G R} C {0} © R © R 
and cap e: R © R —>■ {0} given by 

e = {(u,u,0) I V e R} C R©R©{0}. 

It is straightforward to check these satisfy the zigzag identities. 
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Dagger structure 

Given symplectic vector spaces U, V, observe that the map 

(-)M 7 © y ^ F © c/; 

{u,v) I— 

takes Lagrangian subspaces of the domain to Lagrangian subspaces of the codomain. Thus we can 
view it as a map (—)^ taking morphisms L: U V of LagrRel to morphisms ■. V U. This 
defines a dagger structure on LagrRel, which makes this category into a symmetric monoidal dagger 
category. 

Moreover, every object in LagrRel has a dagger dual: it is clear that r;! = eo a. This category 
thus becomes a dagger compact category. 

6.5 Names of Lagrangian relations 

This brief subsection illustrates a guiding principle of this paper: duals for objects allow us to blur 
the distinction between composition of morphisms and the tensor product of morphisms. We will 
make use of this when we prove the functoriality of the black box functor. 

Observe that a Lagrangian relation L: {0} —?■ R is the same as a Lagrangian subspace of V. 
Moreover, given a Lagrangian relation L: U ^ V, 

U 



V 

compactness allows us to view it as a Lagrangian relation {0} —>■[/© V: 



U L 


We call this subspace the name of the Lagrangian relation L; indeed, we have used this one-to-one 
correspondence between morphisms and their names to define Lagrangian relations. 

By compactness, we have the equation 



Here the right hand side is the name of the composite M o L oi Lagrangian relations, while the left 
hand side is the direct sum of the names of L and M post-composed with the Lagrangian relation 

U V V W 



U W 
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Thus this relation above, the product of a cap and two identity maps, enacts composition of La- 
grangian relations. The proof of Proposition 6.8, that the composite of two Lagrangian relations is 
again a Lagrangian relation, makes use of this fact. We shall also return to it when discussing our 
functor Circ —>■ LagrRel. Note the similarity in form between our diagrams of cospans and of names 
of Lagrangian relations. 


7 Ideal wires and corelations 

In the previous section our exploration of the meaning of circuit diagrams culminated with our 
understanding of behaviors as Lagrangian subspaces. We now turn our attention to how circuit 
components fit together, and the category of operations. In this section we shall see that the algebra 
of connections is described by the concept of corelations, a generalization of the notion of function 
that forgets the directionality from the domain to the codomain. We then observe that Kirchhoff’s 
laws follow directly from interpreting these structures in the category of linear relations. 


7.1 Ideal wires 

To motivate the definition of this category, let us start with a set of input terminals X, and a set of 
output terminals Y. We may connect these terminals with ideal wires of zero impedance, whichever 
way we like—input to input, output to output, input to output—producing something like: 



In doing so, we introduce a notion of equivalence on our terminals, where two terminals are equivalent 
if we, or if electrons, can traverse from one to another via some sequence of wires. Because of this, 
we consider our perfectly-conducting components to be equivalence relations on X + Y, transforming 
the above picture into 



The dotted lines indicate equivalence classes of points, while for reference the grey lines indicate 
ideal wires connecting these points, running through a central hub. 

Given another circuit of this sort, say from sets Y to Z, 
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we may combine these circuits in to a circuit X to Z 






by taking the transitive closure of the two equivalence relations, and then restricting this to an 
equivalence relation on X + Z. This in fact defines a dagger compact category of fundamental 
importance: the category of corelations. 

Ellerman gives a detailed treatment of corelations from a logic viewpoint in [19], while basic 
category theoretic aspects can be found in Lawvere and Rosebrugh [29]. However, to make this 
paper self-contained, we explain them here. 


7.2 The category of corelations 

In the category of sets we hold the fundamental relationship between sets to be that of functions. 
These encode the idea of a deterministic process that takes each element of one set to a unique 
element of the other. For the study of networks this is less appropriate, as the relationship between 
terminals is not an input-output one, but rather one of interconnection. 

In particular, the direction of a function becomes irrelevant, and to describe these interconnec¬ 
tions via the category of sets we must develop an understanding of how to compose functions head to 
head and tail to tail. We have so far used cospans and pushouts to address this. Cospans, however, 
come with an apex, which represents extraneous structure beyond the two sets we wish to specify a 
relationship between. Corelations arise from omitting this information. 

Definition 7.1. A corelation a: X ^ Y between finite sets X and Y is a partition a of the 
disjoint union X + Y. 

That is, for finite sets X and Y, a corelation is a collection of nonempty subsets a = {Ai, A 2 ,..., An} 
of X + Y such that 

(i) a does not contain the empty set. 

(ii) U?=i^^=^ + y- 

(iii) Ai n Aj = 0 whenever i j. 

For example, we can take a circuit of ideal wires with X as the set of inputs and Y as the set of 
outputs: 


X 



V 


and dehne a corelation a: X ^ Y for which terminals lie in the same set of the partition a 
{Hi, A 2 ,..., An} when we can travel from one terminal to another following a path of wires: 



Y 
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The discussion in the previous section then motivates a rule for composing corelations. We 
compose a corelation a: X ^ Y and a corelation Y Z hy finding the finest partition on 
X -\-Y + Z that is coarser than both a and /3 when restricted to X + Y and Y + Z respectively, and 
then restricting this to a partition onX + Z. More explicitly, the corelation jSoa = {Ci,C 2 , ■ ■ ■, Cm} 
is the unique set of pairwise disjoint Ci of the form 

C^=[jA,r^XV^\JBkf^Z 

j k 


for j, k varying over indices such that 


\jA,nY = \jBknY. 
j k 

This rule for composition is associative, as both pairwise methods of composing relations a: X ^ Y, 
p-.Y Z, and y.Z^W amount to finding the finest partition on X + Y + Z + W that is coarser 
than each of a, f3, and 7 when restricted to the relevant subset, and then restricting this partition 
to a partition on X + W. The pictures in the previous section make this clear. 

The identity corelation lx'- X ^ X on a. set X is the map [lx, Ijc ]: X + X ^ X. Equivalently, 
it is the partition ot X + X such that each partition comprises exactly two elements: an element 
X & X considered as an element of both the first and then the second summand oi X + X. We thus 
define a category: 

Definition 7.2. Let Corel be the category with finite sets as objects and corelations between finite 
sets as morphisms. 

This category becomes monoidal under disjoint union of finite sets, using the fact that the union 
of partition ot X + Y and a partition of X' + T' is a partition of (X + X') + (X + Y'). Moreover, 
the commutativity and associativity of the coproduct of finite sets also allows implies the category 
is symmetric monoidal and indeed dagger compact: the braiding and cups and caps are respectively 
given by the partitions of (X + X) + (X + X) and X + X where two elements are in the same part of 
the partition if and only if they are equal as elements of X or X, while a dagger is given by simply 
considering a partition of X + X as a partition of X + X A 

We can get a corelation a: X ^ Y from a surjection /: X + X ^ S' for some set S, by taking 
the sets in the partition a = (Ai,... ,X„) of X + X to be the inverse images of the points of S. 
Two surjections f: X + Y —)-S, / :X + X —>-S' give the same corelation if and only if there is an 
isomorphism g: S ^ S' making the obvious triangle commute: f' = g o f- 

By the universal property of the coproduct, surjections X + Y ^ A are in one-to-one correspon¬ 
dence with jointly epic cospans X —^ X. Thus, corelations can also be seen as isomorphism 
classes of jointly epic cospans. This lets us compose corelations by composing cospans and applying 
a correction. That is: given corelations a: X —>■ X, ,5: X —>■ Z, we may choose cospans representing 
them and compose these cospans. The composite cospan may not be jointly epic. However, we can 
then replace the apex of the composite cospan by the joint image of the feet. The resulting jointly 
epic cospan gives the composite corelation fi o a: X Z. More precisely, we have the following 
proposition. 

Proposition 7.3. There is a strict symmetric monoidal dagger functor 

Cospan(FinSet)-^ Corel 

mapping any cospan to the corelation it defines. 

‘^In fact, if we take a skeleton of Corel we obtain the PROP for ‘extraspecial commutative Frobenius monoids’, as 
defined by Baez and Erbele [4]. We do not need this here, but it helps tie our current work to other work on categories 
in electrical engineering. A proof can be found in Coya and Fong [15]. 
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Proof. This functor maps each finite set to itself. We thus need only discuss how it acts on mor- 
phisms. A cospan in FinSet comprises a pair of functions X ^ N Y. Restricting the apex N 

down to the joint image f{X) U g{Y) gives a jointly epic cospan X ^ f{X) U g{Y) ^ Y, and so a 
corelation X ^ Y. As the elements of the apex not in the image of maps from the feet play only 
a trivial role in the pushout, and hence in composition of cospans, this map is functorial. It is now 
readily observed that the functor is symmetric monoidal. □ 

All this is dual to the more familiar connection between spans and relations, where a relation is 
seen as an isomorphism class of jointly monic spans. This explains the name ‘corelation’. Note that 
neither relations nor corelations are a generalization of the other. The key property of corelations 
here is that it forms a compact category with disjoint union of finite sets as the monoidal product. 
This is not true of the category of relations between finite sets. 

Via the above proposition, the symmetric monoidal dagger functor embedding FinSet into 
Cospan(FinSet) gives rise to a symmetric monoidal dagger functor embedding FinSet into Corel. 
Composing this with the dagger structure on Corel, we also obtain a symmetric monoidal dagger 
functor embedding FinSet°’‘’ into Corel. Corelations thus give a method of composing functions 
regardless of the direction of those functions. Given some not necessarily directed path of functions, 
for example 

A —^ B <— C <— D —> E —> F, 

considering these functions as corelations gives a way to compose them. 

In particular, while this mode of composition is simply composition of functions for two functions 
head to tail, and turning a cospan into a corelation for two functions head to head, for two tail to 
tail functions C ^ D ^ E we compute the composite of cospans 


C 


C 




D 


E 




E 


before restricting the apex to arrive at a surjective function from C + E. This implies the following 
lemma. 


Lemma 7.4. Let 



B 



N 


be a pushout square in FinSet. Then the composites of corelations A ^ P ^ B and A N ^ B 
are equal. 


7.3 Potentials on corelations 

Chasing our interpretation of corelations as ideal wires, our aim for the remainder of this section is 
to build a functor 

S: Corel —> LagrRel 

that expresses this interpretation. We break this functor down into the sum of two parts, according 
to the behaviors of potentials and currents respectively. 
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The consideration of potentials gives a functor <f>: Corel —> LinRel, where LinRel is the symmetric 
monoidal dagger category of finite-dimensional F-vector spaces, linear relations, direct sum, and 
transpose. In particular, this functor expresses Kirchhoff’s voltage law: it requires that if two 
elements are in the same part of the corelation partition—that is, if two nodes are connected by 
ideal wires—then the potential at those two points must be the same. 

This functor is a generalization of the contravariant functor FinSet Vect that maps a set to 
the vector space of F-valued functions on that set. 

Proposition 7.5. Define the funetor 


$: Corel —s- LinRel, 

on objeets by sending a finite set X to the vector space F^, and on morphisms by sending a corelation 
a: X to the linear subspace $( 0 ) o/F^ ©F^ comprising functions f/'y]: X + F —^ F 

that are constant on each element of a. This is a symmetric monoidal dagger functor. 

Proof For coherence maps we take the usual natural isomorphisms F^ ®F^ = F^^^ and {0} = F^. 
We detail only the proof that $ preserves composition; the other properties are straightforward to 
check. 

Let a: X ^ Y and fi: Y Z he corelations. As $ maps corelations to relations, it is enough 
to check both inclusions <I’(/3) o $( 0 ) C $(/? o a) and <I>(/3 o a) C <I>(/3) o $(q;). 

<I>(/3) o $( 0 ) C $(/3 o a): Let </> = [<t>x,4>z] G ‘I’(/3) o <I>(a). We wish to show that for all Ci € fi o a, 
for all c, c' € Ci, we have (j){c) = (j){c'). To this end, note that there exists some fiy'■ F ^ F such 
that (fxY '■= [4>x,4’y] G $(a) and fiyz '■= [4’y,4’z] G Furthermore, by definition this <j)y has 

the property that for all Aj G a, for all a,a' G Aj, we have fixYip) = </>xy(a0) G j3, 

for all b, b' G Bk, we have fiyzib) = fiyzib')- Write fixYZ = [4’x,4’y ^f’z] : A + F + Z F. 

Our desired fact is thus true: for all c, c' G Ci, there exists a sequence c = cq, ci,..., c„ = c' in 
X + Y + Z such that for all £ = 0,1,... n — 1 we have either ci, c^+i G Aj for some j or ci, c^+i G Bk 
for some k, and hence that 

4){c) = (fx+Y+zico) = (j)x+Y+z{ci) = • • • = (j)x+Y+z{Cn) = 4>{c') 


as required. 

$(/? o a) C $(/3) o $(a): Let (j) = [fixyipz] G $(/3 o a). We must show that there exists ^iy : F —> F 
such that [(j)x,f>Y\ G $(a) and [(l)y,(l)z\ G $(/3). We claim 

{ (j){x) if X G A, X, y G Aj for some j; 

(fiyZ) if 0 G Z, y,z G Bk for some k; 

0 if there exist no such x G A or z G Z 

satisfies this. This is well-defined as for all Aj, Aj P X C Ci for some i, so </> is constant on Aj, 
and similarly for all Bj. Thus it does not matter if there are many such x or z with x,y G Aj or 
y,z G Bk- Furthermore, if there exists y such that both x,y G Aj for some x, Aj and y,z G Bk for 
some z, Bk, then the Ci with Aj n A C Ci and Rfc C Z C Ci is unique, so the definitions of fiziy) 
do not conflict. 

Moreover, by construction [<f>x,4>Y] is constant on all Aj, and [4)y,4)z] is constant on all Bk, so 
[(j)x,4'Y] G $(«) and [(/)y,(/)z] G $(/ 3 ) as required. □ 

To recap, we have now constructed a functor <I>: Corel -G LinRel expressing the behavior of 
potentials on corelations interpreted as ideal wires. We now do the same for currents. 
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7.4 Currents on corelations 

Next we consider the case of currents, described by a functor I : Corel —>■ LinRel. This functor now 
expresses Kirchhoff’s current law: it requires that the sum of currents flowing into each part of the 
corelation partition must equal to the sum of currents flowing out. It may also be understood as a 
generalization of the covariant functor Set —>■ Vect that maps a set to the vector space of F-linear 
combinations of elements of that set. 

Proposition 7.6. Define the functor 


I : Corel —5- LinRel 

as follows. On objects send a finite set X to the vector space (F^)*. Given x € X, we write dx 
for the basis vector of (F^)* corresponding to evaluation at x. On morphisms send a corelation 
a: X ^ Y to the linear relation I (a) comprising precisely those 


{ixGy) 


Xydy 

^xeX y^Y 


e (F^)* © (F^)* 


such that for all Ai € a the sum of the coefficients of the elements of Ai n X is equal to that for 
A,nY: 

^ Aa; = ^ Ay. 

x^AiHX y^AiCCY 

This is a symmetric monoidal dagger functor. 

Proof. The coherence maps are the natural isomorphisms (F^)* © (F^)* —>■ (F^+^)* and {0} —> 
(F^)*. Again the only nontrivial task is to check this map I preserves composition. Again we do 
this by checking inclusions I{(3) o 1(a) G I((3 o a) and I{(3 o a) C I(f3) o 1(a). 


1(13) o 1(a) C /(/3 o a): Let (ixAz) G 1(13) o 1(a), with ix = Xxdx and iz = Y^zgz Xzdz. 

Note that this implies that there is some iy = TlyizY Xydy such that (ixAv) € 1(a), (iy, iz) G HP)- 
Then for each Ci G fS o a we have 

^ Aj; = Y, Xx 


xeAjnx 

AjCiXCCi 


yeAjDY 

AjHXGCi 

(By definition of 1(a)) 

yeBkCY 

BknzcCi 

(See composition of corelations) 

Y 

zGBkCiZ 

B^nzcCi 

(By definition of 1(13)) 

Y Z- 

zeCiHZ 



Thus (ix,iz) G 1(13 o a). 
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I{j3 o a) C I[j3) o I{a): The reverse inclusion requires a bit more effort. Let {ixiiz) G /(/? o a). 
We wish to construct some iy = ^y^y ^ydy € (F^)* such that G I{a) and {iy,iz) & liP)- 

ForeachAj G a, write Ifj = nJf Requiring that (zjv, *y) lies in/(a) implies zy G (F^)* 

must satisfy the linear constraints 


Ay — Kj. 

yeAjnv 


In terms of circuits, this says the net current flowing out of the terminals y G Aj fl V is Kj. Note 
that each Ay appears exactly once in the above set of constraints, and so a solution exists. Similarly, 
defining Kj = XzeBjfiz^z each Bj G /3, requiring (zy,zz) G I{j3) implies iy must satisfy the 
^/3 linear constraints 

Ay = Kj, 

yeBjOY 

where again each Ay appears in exactly one of these constraints, and so again a solution exists. 

These two sets of constraints have a simultaneous solution. Indeed, a contradiction would entail 
some subset V oiY such that (i) there exist collections of indices I and /' such that 

F = IJ Ay n r = IJ Bfc n y, 

jei kei' 

and (ii) that summing the Kj over these indices we have 

AjCYCV BkCY<ZV 

Interpreted on circuits, V would be some subset of the intermediate terminals Y on which the net 
current flowing out of a is constrained to differ from the net current flowing into /3. But by the 
definition of composition of corelations, whenever some subset V satifies (i), it is the disjoint union 
of sets of the form Ci fl Y, where Ci G {3 o a —that is, it is the union of connected components of the 
composite circuit. The hypothesis (zx,*y) G I{(3 o a) then implies for each Ci that 

Yi ^3 — y]] — y]] — y]] 

AjnXQCi yeAjHY yeBkHY BkHZCCi 

AjCiXCCi BknZCCi 

SO (ii) cannot hold. This proves the proposition. □ 

Using elementary methods, an algorithm can also be given to construct an explicit solution. 

7.5 The functor from Corel to LagrRel 

We have now defined functors that, when interpreting corelations as connections of ideal wires, 
describe the behaviors of the currents and potentials at the terminals of these wires. In this section, 
we combine these to define a single functor S: Corel LagrRel describing the behavior of both 
currents and potentials as a Lagrangian subspace. 

Proposition 7.7. We define the symplectification functor 

S: Corel —>■ LagrRel 

sending a finite set X to the symplectic vector space 

S{X) =¥^ 
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and a corelation a: X ^ Y to the Lagrangian relation 


S(a) = $(a) © 1(a) C © (F^)* © F^ © (F^)*. 


Then S is a symmetric monoidal dagger functor, with coherence maps inherited from $ and I. 

Proof. As S is the tensor product in LinRel of the symmetric monoidal dagger functors $, /: Corel —> 
LinRel, it is itself a symmetric monoidal dagger functor Corel —>■ LinRel. Thus it only remains to 
be checked that, with respect to the symplectic structure we put on the objects S(X), the image of 
each corelation S(a) is Lagrangian. 

This follows from condition (v) of Proposition 6.4: S(a) is (i) isotropic as, for all (tfxj'ix, 4 ’yPy), 
(4>'x, i'x, f’Y, I-y) S 'S'(c«) we have 


oa{(4>x,ix, 4 >y, *f), (4>'xP'x^ 4>yPy)) 

= ~{d'x{4>x) — ix(4>'x)) + ~ *v(0v) 

= 'i'xif'x) - ixif'Y) + i'xi^'x) - 

= ^ Kdx((l)'x) - ^ Xydyicfy) + ^ Xydyicfy) - ^ X'^^dxicfx) 

xeX y€:Y y& x^X 


= ^ ^ X^dx((j)'x)- ^ Xydy((j)Y) \ + XI ( XI X',^dx((t>x) - X Xydyicfy)] 

AjGa \xeAjnX yeAj-nY j AjGa yeAjCiY j 

= X f X “ X + X f X X 


Aj&a yGA^nV 

= 0 , 


AjGa \a::eAjnX ygA^ny 


and (ii) has dimension equal to 


dim($(a)) + dim(/(a)) = #a + #(X + P) - #a = #(X + Y). 


This proves the proposition. □ 

We have thus shown we do indeed have a functor S: Corel —>■ LagrRel. In the next section 
we shall see that this functor provides the engine of our black box functor, playing the key role 
in showing that we may indeed treat circuit components as black boxes: that is, that circuits that 
behave the same compose the same. Before we get there, we quickly make two relevant observations. 

Example 7.8 (Symplectification of functions). Let f : X —> Y he a function. In this example we 
show that Sf has the form 


Sf = {(f>x,^x,f’Y,^Y) I f>x = rf’Y,^Y = A^x} C e (F^)* eF^ (B (F^)*, 


where /* is the pullback map 


and /* is the pushforward map 



F^; 

<t>° f, 


A-. (F^)* 

*(-) 




49 




(We shall also write /* for the more general map from functions on to functions on that 
takes a function i{—) on F''^ to the function i{— o /).) The claim is then that these pullback and 
pushforward constructions express Kirchhoff’s laws. 

Recall that the corelation corresponding to / partitions X + Y into parts, each of the form 
f~^{y) U {y}. The linear relation $(/) requires that li x € X and y G T lie in the same part of 
the partition, then they have the same potential: that is, 4>x{x) = (pviy)- This is precisely the 
arrangement imposed by f*4>Y = 4>x'- 


(px{x) = (prifix)) = (pviy)- 

On the other hand, the linear relation /(/) requires that if ix = J2 xgx ~ 

then for each y G Y we have 

Aa; = Ay. 

This is precisely what is required by /*: given any (p G F^, we have 


f*ix{(p) 


f* X! ^xdx{cp) = Xxdx{(p o /) = 
x£x x£x y^y 



dy. 


This gives us the above representation oi Sf when / is a function. 

Remark 7.9. We make a remark on our conventions for string diagrams representing the cups and 
caps of the dualities in Corel and LagrRel. 

Note that by the cap duality diagram 


X X 



in Corel we mean the corelation Ux ■ (X + X X 4 — 


0 ), whereas by the cap 


T V 



in LagrRel we mean the Lagrangian relation Uy '■ V (B V —>■ 0 given by the Lagrangian subspace 
{(w,w) & V ®V \v & V}. Although we represent them similarly, the functor S does not map these 
directly onto each other: the image of Ux under S is the Lagrangian relation S'(Ux): F^ © (F^)* © 
F^ © (F^)* ^ 0 given by the Lagrangian subspace 

{{(P,i,(P,-i) I <()GF'^,i G (F^)*} C F-^ © (F-^)* © F^ © (F^)*, 

and in particular a Lagrangian relation F'’'^ © (F^)* © F^ © (F^)* —>■ 0 and not F^ © (F^)* © 
F^ © (F^)* —>• 0. As the cup diagrams in these categories simply denote the dagger image of these 
caps, the analogous statements apply to them too. 

This is an expression of the fact that in cospan categories and Corel there is a canonical self¬ 
duality, while for Lagrangian relations a basis must be picked before there is an isomorphism between 
a symplectic vector space and its conjugate, the conjugate being a canonical dual object. Nonetheless 
as the functor S uses the set X to generate a symplectic vector space, each object in the image of 
S has a canonical symplectomorphism with its dual S*X : F^ © (F^)* —>• F-^ ® (F^)*, with name 
the Lagrangian subspace {(</), i, ip, —i) \ (p G F^, i G (F^)*}.^ The image of the canonical self-duality 

®While we could treat S*X as an atomic notation, we write it here to evoke the idea of the concept of ‘twisted’ 
symplectification introduced in the introduction. 
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Corel is thus given by the composite of this symplectomorphism with canonical duality in LagrRel; 
that is 


( ^ 


S{X) S{X) 


s^x 



Part III 

The Black Box Functor 


We have now developed enough machinery to prove Theorem 1.1: there is a symmetric monoidal 
dagger functor, the black box functor 

■: Circ —LagrRel 

taking passive linear circuits to their behaviors. To recap, we have so far developed two categories: 
Circ, in which morphisms are passive linear circuits, and LagrRel, which captures the external 
behavior of such circuits. We now define a functor that maps each circuit to its behavior, before 
proving it is indeed a symmetric monoidal dagger functor. 


8 Definition of the black box functor 

The role of the functor we construct here is to identify all circuits with the same external behavior, 
making the internal structure of the circuit inaccessible. Circuits treated this way are frequently 
referred to as ‘black boxes’, so we call this functor the black box functor, 

■: Circ —LagrRel. 

In this section we first provide the definition of this functor, and then check that our definition really 
does map a circuit to its behavior. 

8.1 Definition 

It should be no surprise that the black box functor maps a finite set X to the symplectic vector 
space © (F^)* of potentials and currents that may be placed on that set. The challenge is to 
provide a succinct statement of its action on the circuits themselves. To do this, we take advantage 
of four processes we developed in Parts I and 11. 

Let r: X ^ y be a circuit, represented by the decorated cospan 

{X ^Y, r). 

Recall that this means that X and Y are finite sets, T is a F-graph {N, E, s,t,r), and we have a 
cospan of finite sets 


r 



X Y. 


To define the image of T under our functor ■, by definition a Lagrangian relation ■(L) : ■(X) —> 
■ (F), we must specify a Lagrangian subspace 

■ (T) C F^ © (F^)* © F^ © (F’^)*. 
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Recall that to each F+-graph F we can associate a Dirichlet form, the extended power functional 
Pr : F; 

e^E ^ ^ 

and to this Dirichlet form we associate a Lagrangian subspace 

Graph(dPr) = {{<p,d{Pr)<p) | </> S F^^} C F^^ © (F^^)*. 

We consider this Lagrangian subspace as a Lagrangian relation {0} —> F-^ © (F^)*. 

From the legs of the cospan F, the symplectification functor S gives the Lagrangian relation 

o]1': F^ © (F^)* ^ F^ © (F^)* © F^ © (F^)*. 

Writing Y : Y ^ Y for the identity morphism on the finite set Y, S also provides a way of writing 
the identity morphism SY : F^ © (F^)* —>■ F^ © (F^)*. 

Lastly, we have the symplectomorphism 

S*X: F^ © (F^)* — >F^ 

{<j),i) 1—)• ((/), -i). 

The black box functor maps a circuit F to the Lagrangian relation 

{S^X © SY) o S[i, o]1' o Graph(Pr). 

As isomorphisms of cospans of F^-graphs amount to no more than a relabelling of nodes and 
edges, this construction is independent of the cospan chosen as representative of the isomorphism 
class of cospans forming the circuit. 

We picture this as 

■ (X) 


■ (F) = 


■ (V) 

To summarize: 



Definition 8.1. We define the black box functor 

■: Circ —> LagrRel 

on objects by mapping a finite set X to the symplectic vector space 

m{X)=¥^ 

and on morphisms by mapping a circuit T: X ^ Y, represented by the decorated cospan 

{X -^N ■^Y,r) 

to the Lagrangian relation 

■ (F) = {S*X © SY) o S[i, o]1' o Graph(dPr). 
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The coherence maps are given by the natural isomorphisms 


■ (X) 0 ■(y) = © F’^ © (F^)* ^ F^+^ © (F^+’^)* = ■(X + Y) 


and 

{0}^F^©(F^)* = B(0). 

Theorem 8.2. The black box functor is a well-defined symmetric monoidal dagger functor. 

The next, and final, section is devoted to the proof of this theorem. Before we get there, we first 
assure ourselves that we have indeed arrived at the theorem we set out to prove. 

8.2 Minimization via composition of relations 

At this point the reader might voice two concerns: firstly, why does the black box functor refer to 
the extended power functional P and, secondly, since it fails to talk about power minimization, how 
is it the same functor as that defined in Theorem 1.1? These fears are allayed by the remarkable 
trinity of minimization, the symplectification of functions, and Kirchhoff’s laws. 

We have seen that symplectification of functions views the cograph of the function as a picture 
of ideal wires, governed by Kirchhoff’s laws (Example 7.8). We have also seen that Kirchhoff’s laws 
are closely related to the principle of minimium power (Theorems 3.4 and 3.5). The final aspect of 
this relationship is that we may use symplectification of functions to enact minimization. 

Theorem 8.3. Let l : dN N be an injection, and let P be a Dirichlet form on N. Write 
Q = minjv\OAr P for the Dirichlet form on dN given by minimization over N \ ON. Then we have 
an equality of Lagrangian subspaces 

S'i^(Graph((iP)) = Graph(d(5). 

Proof. Recall from Example 7.8 that Si^ is the Lagrangian relation 

= {{cl),L,i,j)OL,i)\(l)€¥^,i e (F^^)*} CF^©(F^)*©F^^©(F^^)*, 

where t*j(0) = ifcfo i), and note that 

Graph(dP) = {(<(), dP^) | <() e F^}. 

This implies that their composite is given by the set 

o Graph(dP) = {{cf o L,i)\(j> ,i € {¥^^)*,dP^ = lJ}. 

We must show this Lagrangian subspace is equal to Graph(dQ)- 

Gonsider the constraint dP^ = L^:i. This states that for all tp G F'^ we have dP^{p) = i{p o (,). 
Letting : X —>■ F be the function sending n G N to 1 and all other elements of X to 0, we see 
that when n G N \ ON we must have 

= dP^iXn) = OXn ol)= i(0) = 0. 

So 4> must be a realizable extension ot tjj = (j> o l. We henceforth write = (j). As /, is injective, 
= (j)o i gives no constraint on if G ¥^^. 

We next observe that we can write Sl^ o Graph(dP) = Graph(dO) for some quadratic form O. 
Recall that Proposition 6.6 states that a Lagrangian subspace L of F^^ © (F^^)* is of the form 
Graph((iO) if and only if L has trivial intersection with {0} © (F^)*. But indeed, if ■0 = 0 then 0 is 
a realizable extension of 0, so = dPo = 0, and hence i = 0. 


dP 

dp(n) 
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It remains to check that O = Q. This is a simple computation: 


0{tp) = dO^{'4i) = oi) = L^dQ^{'ip) = dP^(%i)) = P('(/j) = Q{tp), 

where ip is any realizable extension oi ip £ . □ 

Write t : dN —>■ N for the inclusion of the terminals into the set of nodes of the circuit, and 
: X ON, : y dN for the respective corestrictions of the input and output map to 
dN. Note that [z,o] = to Then we have the equalities of sets, and thus Lagrangian 

relations: 

{S^X © ly) o o]1' o Graph((iPr) = {S*X © SY) o o Sl^ o Graph(dPr) 

= (S'^X © 5'y) o o Graph(dQr) 

U X 

i?€Graph(ciQ) 

We see now that Theorem 8.2 is a restatement of Theorem 1.1 in the introduction. 


9 Proof of functoriality 

To prove that the black box construction is indeed functorial, we factor it into three functors. These 
functors are each symmetric monoidal dagger functors, so the black box functor is too. 

We first make use, twice, of our results on decorated cospans in Section 4, showing the existence of 
categories of cospans decorated by Dirichlet forms and then Lagrangian subspaces, and the existence 
of functors from the category of circuits to each of these. This proceeds by defining functors Dirich 
and Lagr describing decorations and applying Lemma 4.2 to construct the corresponding cospan 
categories, and then by defining the relevant natural transformations and applying Lemma 4.4 to 
construct the desired functors. 

The third functor takes cospans decorated by Lagrangian subspaces, and black boxes them to 
give Lagrangian relations between the feet of such cospans. The functoriality of this process relies 
on interpreting corelations as Lagrangian relations. 

This gives a factorization 

■: Girc-DirichGospan- > LagrGospan->■ LagrRel 

The following subsections deal with these functors in sequence, defining and proving the existence 
of each one using the techniques of Part 11. 

9.1 From circuits to Dirichlet cospans 

To begin, we define a category of cospans of finite sets decorated by Dirichlet forms. This might be 
seen as a solution to our sought after composition rule for Dirichlet forms—but not quite, as it also 
requires the additional data of a cospan. Nonetheless, it has the property that there is a functor 
from the category of circuits to this so-called category of Dirichlet cospans. 

The category of Dirichlet cospans 
Proposition 9.1. Let 

Dirich: (FinSet,©) —> (Set, x) 
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map a finite set X to the set Dirich(X) of Dirichlet forms on X, and map a function f: X ^ Y 
between finite sets to the pushforward function 

Dirich(/): Dirich(X) —)■ Dirich(y); 

Q I—s- f*Q, 

where f^:Q maps (j) S to Q{4>o /)). This defines a functor. 

Moreover, equipping this functor with the family of maps 

5n,m '■ Dirich(A^) x Dirich(M) —Dirich(A^ + M); 

{Qn, Qm) I—>■ Qn{— o ln) + Qm{— o lm), 

where ln ■ N ^ N + M, lm : M ^ N + M are the injections for the coproduct, and with unit 

: 1 —> Dirich(0); 

• I—^ (F®* ^ F; 0 0) 

defines a lax symmetric monoidal functor. 

Proof. Functoriality is just the associativity of composition of functions; lax symmetric monoidality 
is the associativity, unitality, and commutativity of addition of Dirichlet forms. □ 

By Lemma 4.2, we thus obtain a category DirichCospan where a morphism is a cospan of finite 
sets whose apex is equipped with a Dirichlet form. In particular, note that we have overcome 
our inability to define a category whose morphisms are Dirichlet forms. An identity morphism in 
DirichCospan is a circuit whose set of inputs is equal to its set of outputs, with all of its nodes being 
terminals, and no edges: 

• 

inputs 





The functor Circ —DirichCospan 

We have now constructed two symmetric monoidal functors 

(Circuit, p), (Dirich, S): (FinSet, +) —>■ (Set, x) 


that describe the circuit structures and Dirichlet forms we may put on the set respectively. We have 
also seen, motivating our discussion of Dirichlet forms, that from any circuit we can obtain a Dirichlet 
form describing the power usage of that circuit. Now we shall see that this process respects the tensor 
product. More precisely, it specifies a monoidal natural transformation between Circuit and Dirich. 
By Lemma 4.4, this gives us a strict symmetric monoidal dagger functor Circ —>■ DirichCospan. 


Proposition 9.2. Let 

be the collection of functions 

Oat: Circuit(A^) 

{N,E,s,t,r) 


a\ (Circuit, p) 


Dirich(iV); 
[ <() e F^ 


(Dirich, 8) 


2 r(e) 


Then a is a monoidal natural transformation. 
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Proof. Naturality requires that the square 


Circuit (A^) > Dirich(iV) 


Circuit(/) 


Dirich(/) 


Circuit(M) > Dirich(M) 


commutes. Let {N, E, s,t,r) be an F+-graph on N and f:N —>■ M be a function N to M. Then 
both Dirich(/) o aN and ckm o Circuit(/) map (iV, E, s, t, r) to the Dirichlet form 

eGE 

Thus both methods of constructing a power functional on a set of nodes M from a circuit on N and 
a function N ^ M produce the same power functional. 

To show that a is a monoidal natural transformation, we must check that the square 


Circuit(iV) x CircuitDirich(jV) x Dirich(M) 




Circuit (iV + M) ■ 


Ocn + M 


-> Dirich(A^ + M) 


and the triangle 



commute. It is readily observed that both paths around the square lead to taking two graphs and 
summing their corresponding Dirichlet forms, and that the triangle commutes immediately as all 
objects in it are the one element set. □ 


From Lemma 4.4, we thus obtain a strict symmetric monoidal dagger functor 


Q : Circ = CircuitCospan —> DirichCospan. 


Roughly, this says that the process of composition for circuit diagrams is the same as that of 
composition for Dirichlet cospans. Note that although this functor preserves much of the information 
in circuit diagrams, it is not a faithful functor. For example, applying Q to a circuit with edges e, e' 
from the node m to the node n of resistance re and re' respectively, we obtain the same result as 
for the circuit with just one edge e" from m to n whose resistance re" is given by 

1 _ 1 ^ 1 

re" re re' 


9.2 Prom Dirichlet cospans to Lagrangian cospans 

The next step is to show that our process of turning a Dirichlet form into a Lagrangian subspace—by 
taking the graph of its differential—is functorial. 
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The category of Lagrangian cospans 

We begin by describing a category where morphisms are cospans decorated by Lagrangian subspaces. 

Proposition 9.3. Define 

Lagr: (FinSet, +) —(Set, x) 

as follows. For objects let Lagr map a finite set X to the set Lagr(X) of Lagrangian subspaces of 
the symplectic vector space © (F^) . For morphisms, recall that a function f: X ^ Y between 
finite sets may be considered as a corelation, and the symplectification functor S thus maps this 
corelation to some Lagrangian relation Sf: F'’^ © (F'’^’) —>■ F^ © (F^) . As Lagrangian relations 

map Lagrangian subspaces to Lagrangian subspaces (Proposition 6.8), this gives a map: 

Lagr(/): Lagr(X) ^ Lagr(r); 

L^SfiL). 

Moreover, equipping this functor with the family of maps 

Xn,m- Lagr(7V) x Lagr(M) —s- Lagr(iV + M); 

(Ljv, Lm) I—5- © Lm, 


and unit 


Ai: 1 —^ Lagr(0); 

.^{ 0 } 

defines a lax symmetric monoidal functor. 

Proof. The functoriality of this construction follows from the functoriality of S; the lax symmetric 
monoidality from the relevant properties of the direct sum of vector spaces. □ 

We thus obtain a dagger compact category LagrCospan. 


The functor DirichCospan —>• LagrCospan 

We now wish to construct a strict symmetric monoidal dagger functor DirichCospan LagrCospan 
using Lemma 4.4. For this we need a monoidal natural transformation 

(Dirich, (5), (Lagr, A): (FinSet,©) —(Set, x). 


Proposition 9.4. Let 


P: (Dirich, 6) (Lagr, A) 


be the collection of functions 

Pn: Dirich(A^) —> Lagr(A^); 

Q ^ {{f, dQ^) I 0 e F^^} C F^^ © (F^)*. 
Then P is a monoidal natural transformation. 

Proof. Naturality requires that the square 


Dirich(iV) 

Dirich(/) 

Dirich(M) 


pN 


/3m 


■> Lagr(A^) 

Lagr(/) 

■ Lagr(M) 
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commutes for every function f: N —> M. This is primarily a consequence of the fact that the 
differential commutes with pullbacks. As we did in Example 7.8, write /* for the pullback map and 
/* for the pushforward map. Then Dirich(/) maps a Dirichlet form Q on TV to the form f^Q, and 
/3 m in turn maps this to the Lagrangian subspace 

CF^(B(¥^y. 

On the other hand, /Sat maps a Dirichlet form Q on to the Lagrangian subspace 

CF"©(F^)*, 

before Lagr(/) maps this to the Lagrangian subspace 

{{y, UdQy I V- e F^, <^ = ny)) C F^ © (F^y. 


But 

— d(^fj 

so these two processes commute. 

Monoidality requires that the diagrams 


Dirich(A^) x Dirich(M) ^^^^'> Lagr(TV) x Lagr(M) 


Sn,m 
Dirich(iV + M) 


0N + M 


^N,M 

-4 Lagr(TV + M) 


and 


1 



Dirich(0)- 

^0 



Lagr(0) 


commute. These do: the Lagrangian subspace corresponding to the sum of Dirichlet forms is equal 
to the sum of the Lagrangian subspaces that correspond to the summand Dirichlet forms, while 
there is only a unique map 1 —>■ Lagr(0). □ 

From Lemma 4.4, we thus obtain a strict symmetric monoidal dagger functor 

DirichCospan —>• LagrCospan, 

which simply replaces the decoration on each cospan in DirichCospan with the corresponding La¬ 
grangian subspace. 


9.3 Prom cospans to relations 

At this point we have checked that the process of reinterpreting a circuit as a Lagrangian subspace of 
behaviors is functorial. Our task is now to integrate this information as more than just a ‘decoration’ 
on our morphisms. This process constitutes a monoidal dagger functor 

LagrCospan —LagrRel. 

This factor of the black box functor is the one that gives it its name; through this functor we finally 
seal off the inner structure of our circuits, leaving us access only to the behavior at the terminals. 
Its purpose is to take a Lagrangian cospan, which captures information about the behaviors of a 
circuit measured at every node, and restrict it down to a relation detailing the behaviors simply on 
the terminals. 
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Proposition 9.5. We may define a symmetric monoidal dagger functor 


LagrCospan —s- LagrRel 

as follows. On objects let this functor take a finite set X to the symplectic vector space © (F'’^)*. 
On morphisms let it take a Lagrangian cospan 

{X ^ N ^Y-L C¥^ ® (F^)*) 

to the Lagrangian relation 

{S*X © SY) o S[i, o]1' o L C F^©(F^)* © F^ © (F^)*. 

Proof. The coherence maps here are the usual natural isomorphisms 

F^ © (F^)* © F^ © (F'^)* ^ F^+’^ © (F'^+'^)* 


and 

{0} ^F^©(F^)*. 

It is now routine to observe the symmetric monoidality and dagger-preserving nature of this con¬ 
struction, as well as that it preserves identities. Finally, we must check that composition is preserved. 

Using the concept of names introduced in Section 6.5, this comes down to checking this equality 
of Lagrangian subspaces: 



where the left hand side is the composite in LagrRel of the images of the Lagrangian cospans 
{X N Y'L) and (Y M Z'K), and the right hand side is the image of their 
composite in LagrCospan. (Recall that we write Jn ■ N ^ N ©y M and Jm '■ XI ^ N +y M for 
the maps given by the pushout.) 

Recalling Remark 7.9 pertaining to the functor S and duals for objects, this implies that it is 
enough to check the equality of corelations 


N M N M 



X Z X z 


Writing this instead as a commutative diagram, we wish to prove the equality of the two corelations 
N + M ^ X + Y-. 


[3n,3m] 


' X + Y + Y + Z ■ 


.N +Y M 


idx-I-[idy .idy]-l-idz 





^X + Y + Z 
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As the right-hand part admits the factorization 


N -\-Y M - - fTWOix-.lMOOzl 

[iNOix,iNOOY ,3 moozT 



x + z 


it is enough to check the following diagram commutes when interpreted as a pair of morphisms from 
N + M to X + Y + Z in the category of corelations: 

N +Y M 

[jNOix Jnooy Jmooz] 

N+M X+Y+Z 

X + Y + Y + Z 



By Lemma 7.4, this is equivalent to checking that it is a pushout square in FinSet. This is so: the 
square commutes in FinSet as it is the sum along the lower right edges of the three commutative 
squares 



N+y M 

[iiVJML^ jN°OY=jM°iY 


N + M 

OY+iY 


Y + Y 


Y 

''[idy .id-r] 


N +Y M 



and given any other object T and maps /, g such that 



X + Y + Y + Z 


commutes, there is a unique map N +y M ^ T defined by sending a in iV +y M to /(d), where 
d is a preimage of a in iV -f M under the coproduct of pushout maps [jat, Jm]- This is well-defined 
as the preimage of a is either unique or equal to {oy( y), iy (u)} for some element y GY, and the 
commutativity of the above square containing T implies that f{oY{y)) = fi’iviy))- This proves the 
functoriality of the map LagrCospan LagrRel defined above. □ 

The three functors of this section compose to give the black box functor 

■: Circ —)• LagrRel. 

Since they are each separately symmetric monoidal dagger functors, the black box functor is too. 
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